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SECTION  I 
INTRODUCTION 


1.  BACKGROUND 

Fracture  mechanics  presently  provides  the  best  available  tools  to 
quantitati vely  assess  the  influence  of  preexistent  cracks  in  structures. 
Although  fracture  mechanics,  as  we  now  know  it,  is  a  relatively  new  tech¬ 
nology  (since  the  mid  1950s),  the  basic  ideas  were  already  presented  by 
Griffith  (Reference  1)  in  the  early  1920s. 

Current  linear  elastic  fracture  mechanics  (LEFM)  concepts  assume 
that  the  stress  intensity  factor  K,  the  parameter  that  relates  load, 
crack  length,  and  geometry,  controls  fracture  (K  =  =  constant  at 

fracture)  and  crack  propagation  (Reference  23). 

Stress  intensity  factor  calibrations  are  required  for  any  fracture 
or  crack  propagation  analysis.  Since  crack  tip  analysis  can  be  complex, 
experimental  K  solutions  often  become  necessary  to  verify  or  supplement 
analytical  or  numerical  solutions. 

For  complex  type  geometries  and  loading  configurations  which  commonly 
exist  in  aircraft  structures,  experimentation  may  be  the  primary  method 
to  obtain  a  reasonably  accurate  and  dependable  solution.  Conducting 
experiments  on  actual  parts  and  loadings  is  often  a  very  difficult  and 
costly  procedure. 

The  main  motivation  for  this  investigation  was  to  find  a  technique 
that  may,  in  some  cases,  greatly  reduce  the  experimental  effort  required 
in  order  to  obtain  a  K  solution. 

An  analytical  approach  to  determine  the  stress  intensity  factors  K 
was  discussed  by  Bueckner  (Reference  2)  and  Rice  (Reference  3).  They 
showed  that  once  the  displacement  field  and  stress  intensity  factors  are 
known  for  one  geometry  and  loading  configuration  referred  to  subsequently 
as  Case  1,  K  may  be  obtained  for  any  other  loading  (Case  2)  applied  to 
the  same  crack  geometry  (References  4  through  8,  21).  The  method  depends 
mainly  on  the  reciprocal  theorem  and  strain  energy  expressions.  A  load 
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and  material  independent  weight  function  can  be  constructed  to  relate  the 
stress  intensity  factor,  material  properties,  and  crack  length  together 
with  crack  surface  opening  (for  more  details  see  Part  II  and  Appendix  A). 

The  significant  advantage  of  this  procedure  lies  in  the  fact  that 
only  the  first  problem  (Case  1)  needs  to  be  solved  directly.  Construction 
of  the  weight  function  for  this  case  allows  ready  determination  of  K  for 
any  other  loading  on  the  same  crack  geometry  (Cases  2,  3,  etc.).  Stress 
intensity  factors  may  then  be  obtained  for  a  variety  of  other  specimen 
loadings,  some  of  which  may  be  highly  complex  compared  to  the  original. 

In  addition,  the  computational  expenses  are  minimal  because  only  a  single 
relatively  simple  load  configuration  needs  to  be  calculated. 

The  only  additional  information  required  for  any  other  case  is  the 
stress  distribution  at  the  crack  location  for  a  noncracked  body  under  the 
load  configuration  in  question.  This  stress  distribution  is  available, 
in  many  cases,  from  the  static  crack  free  stress  analysis  or  can  be 
obtained  with  much  less  effort  than  a  cracked  body  analysis  requires. 

2.  PURPOSE 

Although  weight  functions  have  been  employed  with  analytical,  finite 
elements  and  other  numerical  methods,  the  author  is  unaware  of  attempts  to 
combine  weight  functions  with  experimental  procedures. 

The  objective  of  this  work  was  to  show  that  the  weight  function  idea 
can  assist  the  experimentalist  by  allowing  him  an  opportunity  to  obtain 
the  stress  intensity  factor  solution  for  a  complex  loading  configuration 
from  the  results  of  a  simple  experiment,  and  also  to  obtain  K  solutions 
for  a  variety  of  loadings  from  a  single  set  of  experiments.  This  will 
greatly  reduce  the  complexity  and  scope  of  the  experimental  effort. 

3.  GENERAL  APPROACH 

The  weight  function  procedure  is  demonstrated  with  experimental 
results  for  an  edge  crack  specimen. 
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The  reference  (Case  1)  problem  has  been  chosen  as  an  edge  cracked 
strip  subjected  to  four-point  bending.  Crack  surface  displacement  measure¬ 
ments  near  the  crack  tip  were  obtained  by  a  laser  interferometric  method 
following  the  technique  developed  at  the  Materials  Laboratory  (References 
9,  10,  19).  These  results  gave  the  mode  I  stress  Intensity  factor  Kj  for 
that  case.  Crack  mouth  opening  measurements,  using  a  clip  gage,  were 
then  used  to  construct  the  crack  surface  displacement  function  incor¬ 
porating  Orange  (Reference  11)  conic  section  approximation.  These 
experimental  results  were  used  to  evaluate  the  weight  function  and  predict 
the  stress  Intensity  factor  Kj  for  other  cases,  which  were  chosen  to  be 
uniform  tension  and  three-point  bending.  The  results  were  compared  with 
the  known  solution  given  in  References  12  through  14. 
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SECTION  II 
THEORY 


1.  GENERAL  ASSUMPTIONS 

a.  Linear  elastic  fracture  mechanics  (LEFM)  applies  i.e.,  the  amount 
of  plasticity  near  the  crack  tip  is  relatively  small 

^<0.1 

where  r^  is  the  radius  of  plastic  zone.  (See  Hertzberg,  Reference  15,  and 
Appendix  B  for  a  more  detailed  discussion.) 

b.  The  problem  is  assumed  to  be  purely  two-dimensional;  i.e.,  no 
variations  through  the  thickness  are  considered. 

c.  Body  forces  are  assumed  to  be  negligible. 

d.  Only  pure  Mode  I  crack  openings  will  be  discussed. 

2.  STRESS  INTENSITY  FACTOR 

The  stress  intensity  factor  K  is  the  LEFM  parameter  that  relates  load, 
crack  length  and  geometry.  The  LEFM  approach  to  predicting  crack  growth 
assumes  that  K  controls: 

a.  Fracture  (K  =  K^.  =  constant  at  fracture) 

b.  Crack  propagation  due  to  fatigue 

c.  Crack  propagation  due  to  stress  corrosion 

Three  modes  of  stress  intensity  factors  are  defined  depending  on  the  crack 
tip  opening  mode  (Figure  1).  They  may  be  stated  as: 


Mode 

I 

Opening  Kj 

=  lim 
r-*-0 

/2-rrr 

0 

y 

(0=0) 

0) 

Mode 

II 

Sliding  Kjj 

=  lim 
r-t-0 

/27tr 

°xy 

(0=0) 

(2) 

Mode 

III 

Tearing  Kjjj 

=  lim 
r-*-0 

/Sir 

cyz 

(0=0) 

(3) 

where  crack  tip  coordinates  are  defined  in  Figure  2. 
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MODE  I 
(Opening) 


MODE  II 
(Sliding) 


MODE  III 
(Tearing) 


Figure  1.  Crack  Opening  Modes 


Figure  2.  Crack  Tip  Stress  Coordinates 
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This  work  deals  only  with  the  opening  Mode  I,  though  it  can  be  easily 
extended  to  other  "modes." 


3.  THE  WEIGHT  FUNCTION 

The  Mode  I  stress  intensity  factor  Kj  under  arbitrary  loading,  shown 
in  Figure  3,  can  be  given  as: 


K 


I 


(4) 


(see  Appendix  A  and  References  2  through  6  for  detailed  discussion), 
where: 


r  is  any  path  chosen  around  the  body  that  includes  all  the 
surface  tractions  and  body  forces, 

A  is  the  region  defined  by  F, 

T  designates  surface  tractions, 

K.  is  the  Mode  I  stress  intensity  solution  corresponding  to  the 
1  loads  T  and  f. 


f  are  body  forces. 


h  =  h(x,y,a) 


_H_  9.n(x»y?a)  weight  function) 


(5) 


n(x,y,a)  is  the  crack  vertical  displacement  for  any  chosen  load 
configuration  (independent  of  T  and  f), 

K,  in  the  weight  function.  Equation  5,  is  the  Mode  I  stress 
1  intensity  solution  that  corresponds  to  the  n  solution,  and 


H 


for  plane  stress 
for  plane  strain. 


The  weight  function  h  is  independent  of  the  loads  T,  f,  and  material 
properties  and  can  be  calculated  for  the  same  geometry  from  another  known 
solution  of  Kj  and  3n/3a. 


i 


Figure  3.  Cracked  Body  Loading  Configuration 


Using  the  principle  of  superposition  of  linear  elasticity,  the 
calculations  of  the  stress  intensity  factor  considering  the  case  of  a 
cracked  body  under  a  certain  loading  configuration  is  equal  to  the  case 
of  the  cracked  body  where  only  the  crack  surface  is  subjected  to  tractions 
p(x)  (coordinates  are  chosen  so  that  the  crack  is  in  x  direction).  Here 
p ( x)  is  the  stress  ay  that  will  occur  at  the  crack  location  for  an 
uncracked  body  subjected  to  the  original  loading  configuration.  This  idea 
is  demonstrated  in  Figure  4. 


If  one  chooses  the  path  T  as  any  path  that  includes  the  crack  surface 
and,  for  simplicity,  takes  a  straight  crack  in  the  x-direction  as  in 
Figure  4  (the  last  assumption  is  not  essential  and  is  used  only  to  simplify 
the  equations),  then  Equation  4  for  Kj  can  be  stated  as: 


K 


I 


h(x,a)p(x)dx 


h(x,a)p(x)d>: 


(6) 


Reversing  the  integration  limits  and  adding  the  two  integrals  we  get: 


K 


I 


h(x,a)p(x)dx 


(7) 
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Figure  4.  Superposition 


K 


I 


+  K, 


© 


Part  B  can  be  viewed  as  a  noncracked  body;  therefore , 


K 


I 


<D 


0 


© 


The  unknown  stress  intensity  solution  for  the  loading  p(x)  is  con¬ 
sidered  Case  2,  and  the  known  solution,  from  which  the  weight  ^unction 
will  be  computed,  is  taken  as  Case  1.  Then  Equation  7  can  be  written  as 


and  from  Equation  5 


KI@=  2/ 


hQfx,a)  = 


(8) 


(9) 
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If  Equation  9  is  substituted  into  Equation  8,  the  result  is: 


(10) 


where 


K_  is  the  Mode  I  stress  intensity  factor  to  be 

A©  found, 

P(x)0  is  the  stress  at  the  crack  location  for  a 

noncracked  body  with  loading  configuration  of 
Case  2, 


K 


is  the  Mode  I  stress  intensity  factor  for 
Case  1  and  will  be  determined  experimentally, 
and 


3n 


is  the  derivative  with  respect  to  crack  length 
of  the  crack  opening  profile  for  Case  1.  n  is 
also  determined  experimentally. 


4.  K  CALIBRATION  BY  CRACK  TIP  DISPLACEMENT 

The  crack  tip  displacement  field  is  defined  by  the  following  relations 
that  represent  the  elasticity  solution  for  the  stress-displacement  field 
in  terms  of  K  near  the  crack  tip  of  a  cracked  body.  One  can  take  these 
displacement  expressions  for  Mode  I  opening  (References  17,  18,  10)  as: 


u. 


u 


K_  ( -  .  0  T  2  j  0 1 

■  sin  *  b*  ■ cos  *1 


(ID 

(12) 
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for  plane  stress;  and 


for  plane  strain. 

Consider  the  vertical  displacement  uy  along  the  crack  surface  for 

which 

0  =  IT 
r  =  t 


where  t  is  a  coordinate  along  the  crack,  starting  at  the  crack  tip.  One 
should  note  to  conform  with  usual  practice,  n  replaces  the  displacement 
in  y  direction  uy.  Recall  from  elasticity  that: 


G  2 ( 1+v) 


(15) 


Substituting  into  Equations  12  and  14  the  crack  displacement  for  plane 
stress  and  plane  strain  can  be  written  as 


The  equation  holds  only  very  near  the  crack  tip 

t  «  a 

and  allows  one  to  calibrate  Kj  by  measuring  crack  opening  at  a  known 
distance  very  near  the  crack  tip. 
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5.  CRACK  OPENING  CALIBRATION 

In  order  to  find  the  slope  with  respect  to  the  crack  length  3n/3a 
of  the  crack  opening  profile,  an  approximation  has  been  used.  Orange 
(Reference  11)  suggested  a  conic  section  approximation  for  a  finite 
width  edge  cracked  plate.  This  approximation  is  based  on  assuming  a 
general  conic  section  function  for  the  crack  opening  shape  and  fitting 
the  crack  tip  radius  of  curvature  and  crack  mouth  opening  to  determine 
the  unknown  coefficients.  Further  development  of  this  relationship  is 
presented  in  Appendix  E.  Following  Orange  and  substituting  for  the 
crack  tip  curvature  in  terms  of  the  stress  intensity  factor,  using 
crack  tip  displacement  equations,  the  following  expression  is  derived: 


(18) 


where  nQ  is  the  crack  mouth  displacement  and  m  is  the  conic  section 
coefficient. 


Y  is  a  nondimensional  stress  intensity  factor 

»  _  KI  (20) 

o/a~ 

Crack  displacements  of  finite  single  edged  cracked  strips  under 
bending  or  tension  generated  by  Orange  (Equation  18)  were  compared  to  the 
collocation  method  presented  by  Gross  (Reference  20)  and  found  to  match 
within  3.2  percent.  The  Orange  equation  was  compared  to  the  Rice 
(Reference  3)  solution  for  a  center  crack  considering  an  infinite  plate. 
The  results  were  found  to  be  exactly  the  same.  The  Orange  equation  was 
also  used  successfully  by  Grandt  (References  6,  21)  for  cracks  in  fastener 
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holes  and  for  radially  cracked  rings.  Substituting  the  expression  for 
Y  Equation  20,  and  m  Equation  19  into  Equation  18,  we  obtain 


(-t)' SJW 

"Usd 


r8Kilt  ♦ 

iv 

eV] 

-  irH2  J 

a2 

irH2a  J 

In  order  to  simplify  the  above  equation. 


Thus, 


V  8KI2 

a2  irH2a 


T)  =  [Bt  +  Rt2] 


Also, 


|J  -  ^  [Bt  +  RT2]  **  *  Jg  iBt  +  Rt2] 
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Figure  5.  Edge  Crack 


3t 


From  Figure  5,  if  t  =  a-x,  then  rr  =  1 . 


Thus, 


U  -  5  IBt  +  Rt’l  iff  t  +  B  +  |S  t!  +  2Rt] 


(27 


where : 


3B 

3a 


16K][ 

ttH2 


(28 


and 


3R 

3a 


ttH2  a2 


16Kj 

irH2a 


3K, 


3a 


(29 


Equation  27  therefore  defines  3n/3a  in  terms  of  the  crack  mouth  displace¬ 
ment  n  and  the  stress  intensity  factor  KT  and  their  derivative  with 
o  i 

respect  to  the  crack  length. 

6.  Kj  FOR  CASE  2 

Substituting  Equation  27  into  Equation  10  and  also  changing  the 
integration  variable  from  x  to  t  one  obtains 

(t^  [Bt+Rt2]  [B+(||+ 2R)t  +  ||t2]dt 
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3  D 

where  B,  R*  ■  and  are  defined  previously  and  calculated  for  Case  1. 
aa  da 

From  Equation  30  it  is  clear  that  measuring  the  stress  intensity 
factor  Kj  and  the  crack  mouth  opening  q  for  various  crack  lengths  and 
calculating  their  derivatives  with  respect  to  "a"  allows  one  to 
construct  the  weight  function  for  the  Case  1  loading  configuration.  By 
knowing  the  stress  distribution  at  the  crack  location  for  a  noncracked 
body  under  the  loading  configuration  of  Case  2  p(t)^  ,  one  may  obtain 
Kj  for  Case  2. 

Computing  this  integral  is  quite  tedious  though  possible  for  very 
simple  forms  of  p(t).  However,  numerical  integration  is  possible  for 
any  p(t)  even  if  represented  by  a  discrete  numerical  solution.  In  per¬ 
forming  numerical  integration,  one  should  be  cautioned  to  the  fact  that 
the  weight  function  and  therefore  the  whole  expression  under  the  integral 
is  square  root  singular  at  the  crack  tip  t=0  but  still  it  can  be  shown 
that  for  the  problems  studied  here  a  limit  exists  for  the  integral  and, 
in  most  cases,  the  numerical  integration  converges  quite  rapidly. 
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SECTION  III 

EXPERIMENTAL  TECHNIQUE 

1.  THE  LASER  INTERFEROMETRY  METHOD 

Stress  intensity  factor  K  calibrations  obtained  by  measuring  crack 
tip  displacements  require  a  technique  capable  of  accurately  measuring 
extremely  small  displacements  very  near  the  crack  tip.  As  a  rule  of 
thumb,  the  relations  for  the  crack  tip  displacement  field  as  expressed 
in  Equations  11  to  14  are  regarded  to  be  sufficiently  accurate  within  a 
distance  of  a/20  or  smaller  from  the  crack  tip,  where  a  is  the  crack 
length.  The  laser  interferometry  technique  described  here  was  reported 
in  References  9,  10,  and  19  to  measure  displacements  of  0.1  microns  at  a 
distance  of  50  microns  from  the  crack  tip. 

Results  were  reported  leading  to  Kj  values  that  varied  up  to  15 
percent  from  the  theoretical  value  (for  individual  points).  Most  of 
these  deviations  were  attributed  to  the  imperfect  crack  and  material 
characteristics  and  not  to  measurement  errors.  The  technique  is  basically 
similar  to  conventional  interferometry  but  only  measures  in-plane 
displacements  on  the  surface. 

Two  small  indentations  are  placed  with  a  diamond  indentor  on  both 
sides  of  the  crack  near  the  crack  tip.  Those  indentations  are  typically 
square  based  pyramids  with  a  base  length  of  20  to  40  microns  and  located 
50  microns  apart  on  both  sides  of  the  crack  (Figure  6).  A  laser  coherent 
source  impinges  upon  the  indentations.  The  beam  is  diffracted  back  at  an 
angle  a  with  respect  to  the  incident  beam  as  shown  schematical ly  in 
Figure  7. 

Since  the  indentations  are  placed  close  together,  the  respective 
diffraction  beam  overlaps  resulting  in  interference  fringe  patterns  on 
either  side  of  the  incident  laser  beam.  A  photograph  of  a  typcial  inter¬ 
ference  pattern  is  shown  in  Figure  8. 
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Figure  8.  Typical  Interference  Fringe  Pattern 

The  relation  between  the  indentation  spacing  and  the  fringe  order 
is  set  by  the  following  equation  (Reference  9): 

d  sin  aQ  -  mX  (31 ) 

where  d  is  the  distance  between  indentations,  X  is  the  wave  length  of 
the  source,  and  m  is  the  fringe  order. 

As  a  load  is  applied  to  the  specimen,  the  crack  opens  and  the  dis¬ 
tance  between  the  indentation  changes.  The  relation  between  the  change 
in  d  and  the  change  in  the  fringe  order  6m  from  a  fixed  point  of  view  is 

fid  =  5m  X/sin  a 

o 

Now  observing  the  fringe  pattern  from  a  fixed  point  allows  one  to  count 
the  number  of  fringes  passing  and  hence  determine  the  change  in  the 
fringe  order.  Equation  32  will  then  give  the  change  in  the  distance  d 
between  the  indentations  that  correspond  to  the  crack  opening.  Averaging 
the  left  and  right  fringe  number  will  eliminate  free  body  motion. 
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2.  TEST  APPARATUS 

The  test  specimens  used  were  flat  strips  of  Aluminum  7075  T651. 
A  small  V  notch  was  placed  on  the  edge  of  the  specimen  to  allow  pre¬ 
cracking  and  a  shallow  groove  was  cut  at  the  edge  as  to  fit  the  clip 
gage  used  to  measure  crack  mouth  openings.  (Figure  9). 


Figure  9.  Test  Specimen 


The  specimen  was  placed  in  a  four-point  bending  fixture.  Load  was  applied 
by  a  standard  Instron  machine  equipped  with  a  compression  load  cell 
(Figure  10). 


Bending  fixture 


Specimen 


Clip  gage 


Load  cell 


Figure  10.  Specimen  Set  Up  in  Bending  Fixture 
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The  crack  mouth  opening  was  measured  by  a  clip  gage  placed  In  the 
premachined  groove.  The  clip  gage  consists  of  two  small  steel  plates  on 
which  strain  gages  were  bonded  (Figure  11).  Precallbrating  the  clip  gage 
allows  one  to  relate  the  mouth  opening  with  the  changes  In  strain  measured. 

Two  pyramid-type  Indentations  have  been  placed  near  the  crack  tip 
by  a  standard  Lietz  microhardener  tester  with  a  diamond  Indentor  of  a 
square  based  pyramid  shape  with  face  angles  of  136  degrees.  The  indentations 
were  impinged  by  a  Spectra  Physics  Model  120  5  MW  HeNe  laser.  The  laser 
beam  wave  length  X  was  0.6328  microns  and  had  a  divergence  angle  of  0.71 
mr  and  a  beam  diameter  of  0.57  mm.  The  laser  was  located  so  that  the 
beam  impinged  perpendicular  to  the  specimen  surface.  The  angle  of  the 
diffracted  beam  aQ  can  be  calculated  directly  by  knowing  the  face  angle 
of  the  indentation.  (Figure  11). 


Figure  11.  Clip  Gage  Setup 

Li 


AFWAL-TR-80-4001 


The  fringes  created  by  the  interference  of  the  reflected  beams  were 
picked  up  by  two  photoresistors  that  were  located  in  a  fixed  place  on  the 
fringe  pattern. 

The  active  face  of  the  photo  resistors  was  masked  leaving  a  narrow 
slit  (smaller  than  the  fringe  spacing)  to  allow  effective  distinction 
between  individual  fringes.  A  photograph  of  the  test  apparatus  is  seen 
in  Figure  13. 

As  load  was  applied  and  the  fringes  started  moving,  the  photoresistor 
created  an  electrical  signal  that  after  proper  amplification  could  be 
recorded  with  the  aid  of  a  strip  chart  recorder. 

3.  TEST  PROCEDURE 

Two  Aluminum  7075  T651  specimens  were  each  used  to  test  5-10  different 
crack  lengths.  The  first  specimen  was  used  for  longer  cracks  and  the 
second  for  shorter  cracks  with  some  overlapping  between  the  regions  that 
allowed  comparing  between  the  two.  The  specimens  were  precracked  in 
the  long  transverse  (LT)  direction  of  the  material  in  three-point  bending 
on  a  Schenck  fatigue  machine  using  cyclic  loads  at  a  frequency  of  20  Hz. 

The  precracking  stress  intensity  factor  varied  from  6000-9000  psi  /Tn . 

The  higher  Kj  was  used  to  initiate  and  propagate  the  crack  to  about  80 
percent  of  the  desired  length,  and  the  lower  Kj  levels  were  applied  to 
get  the  final  sharp  crack. 

Two  indentations  were  placed  on  the  surface  50-200  microns  behind 
the  crack  tip  as  explained  previously.  A  load  of  200  grams  on  the 
indentor  diamond  head  was  found  most  suitable  to  create  good  reflecting 
indentations  with  a  base  square  dimension  of  approximately  30  microns. 

The  specimen  was  mounted  on  the  loading  fixture,  and  the  two  photo¬ 
resistor  output  signals  together  with  the  clip  gage  and  load  cell  signals 
were  monitored  by  a  four-channel  strip  chart  recorder  to  create  a  common 
time  basis. 

The  specimens  were  loaded  typically  to  values  of  Kj  =  12-14  ks 1  /Tn 
so  the  amount  of  plasticity  in  the  specimen  was  kept  small  (Appendix  B). 

Two  load  cycles  were  applied  prior  to  actual  testing  so  that  the  crack  would 
appear  to  its  full  length  on  the  surface  and  would  easily  be  seen  under  a 
microscope. 
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During  the  preload  cycle  the  photoresistors  were  relocated  to  a 
position  such  that  the  maximum  fringe  intensity  could  be  observed.  Since 
the  photoresistors  were  fixed  in  position  while  the  specimen  with  the 
indentation  was  displacing  due  to  the  applied  load,  it  became  necessary 
to  establish  the  final  location  of  the  resistors  under  a  mean  load  so 
that  the  indentations  could  still  stay  within  the  beam  boundaries  over  the 
extreme  deflection. 

Three  load  cycles  were  applied  to  each  specimen  for  each  crack 
length.  After  each  test  the  crack  length  was  measured  in  order  to  account 
for  any  crack  propagation  during  the  test  cycle.  Then  the  crack  was 
propagated  to  the  next  stage  in  the  Schenck  fatigue  machine.  In  addition, 
for  two  different  crack  lengths  the  crack  surface  displacement  was 
measured  along  the  crack  at  2-6  different  locations  using  previously 
applied  indentations  in  order  to  verify  the  Orange  (Reference  11)  conic 
section  approximation  for  the  crack  opening. 
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SECTION  IV 

DATA  REDUCTION  AND  NUMERICAL  TECHNIQUE 
1.  Kj  FOR  CASE  1 

The  fringe  patterns,  crack  mouth  opening,  and  loads  were  all  recorded 
on  a  four-channel  strip  chart  recorder.  A  typical  chart  for  a  single 
loading  cycle  is  shown  in  Appendix  C,  Figure  C-l. 

As  explained  previously,  two  sets  of  fringes  were  created  by  the 
reflected  beam  from  the  indentation.  For  each  set  of  fringes,  a  slope 
of  number  of  fringes  per  unit  load  6m/p  can  be  obtained.  These  slopes  were 
averaged  to  eliminate  free  body  motion.  A  more  detailed  description  of 
the  laser  interferometer  data  reduction  procedure  can  be  found  in 
Appendix  C. 

Using  Equation  32 


n  =  1  6d  =  1  =  1  6m  .  X 
P  2  P  P  2  P  sin  a 

o 

where  X/sin  aQ  is  the  calibration  factor  of  the  interferometer  used  here. 
Equation  33  allows  one  to  compute  the  displacement  per  unit  load.  The 
loading  and  unloading  results  of  three  loading  cycles  were  averaged  to 
give  an  average  displacement  per  unit  load.  The  distance  t  of  the 
indentation  behind  the  crack  tip  was  computed  by  averaging  the  initial 
distance  before  the  test  with  the  one  measured  after  three  test  cycles. 
Since  only  three  cycles  were  applied  and  the  amount  of  crack  propagation 
was  very  small  (a  few  microns),  linear  interpolation  was  adequate  to  give 
a  sufficiently  accurate  estimate. 

Kj  for  Case  1  can  now  be  determined  from  the  crack  displacement  near 
the  crack  tip  (Equations  12  or  14).  In  most  cases  G  is  close  enough  to 
180°  so  the  simplified  form  of  Equation  17  can  be  used. 

Since  the  measurements  are  on  the  surface,  the  stress  field  should 
be  considered  as  plane  stress  and  H  should  be  given  the  value  of  E 
(lxl 0^  for  aluminum  7075T651).  Now  Equation  17  takes  the  form 
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P 


(34) 


This  data  will  usually  be  converted  to  some  nondlmenslonal  form. 

2.  CRACK  MOUTH  OPENING 

The  crack  mouth  opening  nQ  was  recorded  directly  from  the  clip  gage 
readings.  The  crack  mouth  opening  per  unit  load  was  then  readily  obtained 
from  a  nQ  versus  P  plot  using  only  the  linear  part  of  the  curve  and 
neglecting  the  initial  nonlinear  part  which  is  associated  with  crack 
closure  effects. 

3.  NUMERICAL  TECHNIQUE  FOR  Kj  CASE  2  COMPUTATION 

The  rest  of  the  computation,  which  involved  calculating  stress 
intensity  factors  for  different  Case  2  loadings  were  programmed  and 

handled  numerically.  The  KT  data  was  fitted  by  a  suitable  interpolation 

*1® 

curve.  Interpolated  values  of  K.  _  and  — 5^-^  were  computed  for  any 

1q  3a 

desired  crack  length  a/w.  In  a  similar  way,  the  crack  mouth  opening 
data  for  Case  1  was  fitted  by  an  interpolating  curve  that  allows 

one  to  compute  values  of  n  and  for  the  desired  points.  These  inter- 
polations  are  calculated  in  a  different  subroutine  so  any  desired  inter¬ 
polation  function  can  be  tried  without  changing  the  main  program.  A 
comparison  table  for  Kj  and  nQ  for  interpolated  data  versus  actual  data 
is  printed  out  so  one  can  judge  the  accuracy  of  the  interpolation  process 
chosen.  Fourth  order  least  square  polynomials  will  give  reasonable 
results  in  many  cases. 

Using  Orange  (Reference  11)  conic  equation  discussed  in  Appendix  E 
for  the  crack  opening  profile,  the  derivative  of  the  crack  opening  along 
the  crack  with  respect  to  the  crack  length  can  be  computed  through 
Equation  27. 

Since  different  factors  are  used  to  nondlmenslonal Ize  the  K,  and  n 

Kt  I  0 

data,  care  should  be  taken  to  bring  all  the  K. ,  3^  ,  n  and 

3no  ° 

-r-  to  a  common  basis.  Using  Equation  9,  the  weight  function  is 

nd 
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computed  for  each  crack  length  at  any  point  along  the  crack.  The  constant 
H  is  determined  depending  on  plane  stress  or  plane  strain  conditions 
(Appendix  B).  The  stress  distribution  along  the  crack  location  for  a 
noncracked  body  p(t)  is  supplied  by  a  separate  subroutine  and  can  be 
easily  fitted  with  available  data  or  some  known  analytical  solution.  Kj 
for  Case  2  now  is  computed  through  integrating  Equation  30. 

The  integral  is  performed  numerically  using  a  Romberg  integration 
scheme  (Reference  24).  Although  the  weight  function  and  thus  the 
integrand  in  Equation  30  is  singular  at  the  crack  tip,  the  integral  will 
still  converge  to  a  distinct  value  quite  rapidly  for  most  p(t)  functions. 
Because  of  the  crack  tip  singularity,  the  lower  limit  of  this  integral 
cannot  be  set  to  zero.  A  value  close  to  zero  should  be  taken  and  this 
value  can  be  made  smaller  and  smaller  by  an  iterative  process  until  the 
value  of  the  integral  converges  within  a  specified  accuracy. 

Since  the  numerical  integration  process  may  be  executed  many  times 
until  convergence  is  obtained,  an  efficient  algorithm  should  be  used  to 
reduce  the  required  computer  time.  Romberg  integration  scheme  used  here 
is  a  powerful  and  efficient  numerical  integration  technique.  It  is 
based  upon  the  trapezoidal  rule  combined  with  Richardson  extrapolation. 
More  details  of  this  scheme  can  be  found  in  Hornbeck  (Reference  24). 

Using  this  scheme  reduced  the  computer  time  by  a  factor  of  about  200 
compared  to  the  Simpson’s  rule  Integration  scheme.  A  detailed  description 
of  the  computer  program  is  presented  in  Appendix  D. 
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SECTION  V 

EXPERIMENTAL  RESULTS 
1.  EXPERIMENTAL  MEASUREMENTS 

Two  Aluminum  7075  T651  edge  cracked  strips  were  used  in  the  experi¬ 
ments.  The  specimens  were  cut  out  from  a  rolled  plate  and  the  crack  was 
grown  in  the  long  transverse  (LT)  direction  of  the  material. 

The  dimensions  of  the  specimens  and  the  region  of  crack  length  in 
which  they  were  tested  are  shown  in  Figure  14. 


*N  dimension  determined  for  clip  gage  purposes 


Figure  14.  Specimen  Dimensions 
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Specimen  A  was  tested  with  ten  different  crack  lengths  with  crack 
length  to  width  ratios  ranging  from  .2328-. 7517.  Specimen  B  was  tested 
with  five  different  crack  lengths  where  a/w  ranged  from  .1034-. 476. 

Both  specimens  were  loaded  in  a  four-point  bending  fixture  as  shown 
in  Figure  15. 


£ 


Figure  15.  Loading  Fixture  Dimension 


Each  specimen  was  loaded  through  three  loading  cycles,  and  near  crack  tip 
displacement  (by  laser  interferometry)  and  crack  mouth  openings  were 
recorded.  Table  1  represents  the  average  results  for  all  crack  lengths 
measured. 


2.  DATA  INTERPOLATION 

The  stress  intensity  data  calculated  from  the  interferometry 
measurements  was  converted  to  a  nondimensional  form  of 

Kj  KjbW2 

o/ifa  6M/ira 

P  1" 

where  M  is  the  bending  moment  M  =  . 


(35) 


The  crack  mouth  opening  data  was  converted  to  a  nondimensional  form 


n  bW2H 
o _ 

12Ma 


(36) 
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TABLE  1 

EXPERIMENTAL  MEASUREMENT  RESULTS 


Specimen 

Nondimensional 
crack  length 

Distance  of 
indentation 
from  crack 
tip 

No  of  fringes 
per  unit  load 

Crack  mouth 
displacement 
per  unit  load 

a/w 

r 

6m/p 

Vp 

Microns 

1/U>  x  10-3 

in/£b  xlO-6 

B 

.1034 

101 

5.093 

0.568 

B 

.1586 

83.5 

6.623 

0.7453 

B 

.  2064 

88 

7.6713 

0.9836 

A 

.2328 

121.5 

7.3128 

0.856 

B 

.  2646 

86.5 

9.116 

1.268 

A 

.2666 

99.5 

j 

7.580 

1.14 

A 

.3094 

143 

10.180 

1.405 

A 

.346 

56 

7.340 

1.705 

A 

.3978 

72 

9.140 

2.358 

A 

.436 

46.5 

8.0146 

2.727 

B 

.476 

47 

12.808 

4.121 

A 

.4893 

50 

10.130 

3.6585 

A 

.5780 

111.5 

18.354 

6.04 

A 

.6467 

74 

20.460 

9.216 

A 

.7517 

70 

35.500 

20.681 
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The  nondimensional  data  for  Kj  and  nQ  was  interpolated  by  a  least  square 
fourth  order  polynomial.  Table  2  and  Table  3  present  the  actual  experi¬ 
mental  and  the  interpolated  data  for  Kj  (nondimensional)  and  nQ 
(nondimensional)  respectively  including  the  relative  difference  in  percent 
of  the  interpolating  polynomial  at  the  test  point. 

Figure  16  presents  the  actual  experimental  nondimensional  Kj  com¬ 
pared  to  the  known  reference  results  from  Tada  (Reference  12).  Figure  17 
presents  the  actual  experimental  nondimensional  data  for  nQ  compared  to 
the  known  reference  results  from  (Reference  12).  The  average  difference 
between  the  reference  values  of  Kj  and  the  experimental  data  is  6%  while 
the  maximum  difference  for  a  single  measurement  is  up  to  14%.  These 
results  are  similar  in  accuracy  to  the  ones  obtained  by  Macha  (Reference 
10)  for  the  laser  interferometry  technique.  The  values  of  Kj  from  the 
interpolating  polynomial  that  were  actually  used  to  compute  the  weight 
function  are  similar  in  accuracy  with  a  maximum  difference  of  13%  and  an 

average  difference  of  6%  compared  to  Reference  12.  One  should  note 

a  K  t 

that  the  relatively  large  difference  in  the  slope  of  -gj  for  some  of 
the  short  crack  lengths  (Figure  16)  will  have  only  a  minor  influence  on 
the  weight  function  computation  as  can  be  seen  from  examining  Equation  30. 
nQ  data  as  well  as  the  interpolating  polynomial  fits  the  Reference  12 
values  much  better  as  seen  in  Figure  17. 

3.  CRACK  PROFILE  MEASUREMENTS 

For  two  different  crack  lengths  on  specimen  A  the  crack  opening  was 
measured  along  the  crack  at  several  points  in  addition  to  the  point  near 
the  crack  tip  and  at  the  crack  mouth.  These  measurements  were  done  by  the 
laser  interferometry  technique.  The  existing  indentations  left  behind 
as  the  crack  propagated  were  used  for  these  measurements.  The  crack 
openings  along  the  crack  were  compared  to  the  Orange  (Reference  11)  conic 
sections  and  found  to  match  this  approximation  excellently  (Figures  18 
and  19). 
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TABLE  2 

Kj  (NONDIMENSIONAL)  COMPARISON  OF  ACTUAL  DATA  TO 
INTERPOLATED  DATA 


Actual  Kj. 

Interpolated 

Relative 

Difference 

a/W 

KjbW2 

K_bW2 
. 1  * 

in  % 

aM/rra 

aM/ifa 

.  1034 

.9073 

.9467 

4.34 

.1586 

1.0575 

.  9776 

-7.  56 

.2064 

1.0572 

1.0545 

-.26 

.  2328 

1.1137 

1.1060 

-.69 

.2646 

1.1007 

1. 1705 

6.35 

.2666 

1.1862 

1.1746 

-.98 

.3096 

1.2366 

1.2590 

1.81 

.  3460 

1. 3490 

1.3244 

-1.82 

.3978 

1.3804 

1.4086 

2.04 

.4360 

1.4388 

1.4696 

2.14 

.  4760 

1.5646 

1.5424 

-1.42 

.4893 

1.6551 

1.5710 

-5.08 

.5780 

1.8481 

1.8774 

1.59 

.6467 

2.3460 

2.3688 

.  97 

.7517 

3.9513 

3.9413 

-.25 

♦The  interpolating  polynomial  is: 

K  bW2 

— -  =  70. 75 (a/W)  4  -92.45(a/W)  3+42. 95 (a/W) 2-6.53(a/W)+1.257 

aM/rfa 
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TABLE  3 

nQ  (NONDIMENSIONAL)  COMPARISON  OF  ACTUAL  DATA  TO 
INTERPOLATED  DATA 


Actual  n 

o 

Interpolated  nQ 

Relative 
Difference 
in  % 

a/W 

n  bWzH 
o 

"obw!«  . 

12Ma 

12Ma 

.1034 

1.9440 

1.9500 

.  31 

.1586 

1.6643 

1.6491 

-.91 

.  2064 

1.6855 

1.6334 

-3.09 

.2328 

1.6163 

1.6879 

4.43 

.2646 

1.6950 

1.7933 

5.80 

.2666 

1.8796 

1.8011 

-4.18 

.  3096 

1.9961 

1.9923 

-.19 

.3460 

2.1660 

2.1823 

.75 

.3978 

2.6055 

2.4945 

-4.26 

.4360 

2.7493 

2. 7696 

.74 

.4760 

3.0632 

3.1251 

2.02 

.4893 

3.2866 

3.2648 

-.66 

.  5780 

4.5933 

4.6431 

1.08 

.6467 

6.6244 

6.5725 

-.78 

.7517 

12.0980 

12.1082 

.08 

*The  interpolating  polynomial  is: 


n  bWH 

=  210.38  (a/W)  l'-282.57(a/W)  3  +  150.74  (a/W)  2-32.16(a/W)  + 


Tada  (Reference  12) 


—  Interpolation  curve 
for  the  experimental 
data 

^  Experimental  data. 
Specimen  A 

A  Experimental  data, 
Specimen  B 
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Orange  Equation  (Reference  11) 


Figure  18.  Nondimensional  Displacement  Along  Crack  Surface  Compared 
to  Orange  Equation  (Reference  11)  Pure  Bending  a/W=0.436 


to  Orange  Equation  (Reference  11)  Pure  Bending,  a/W=0.7517 
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4.  CASE  2  LOADINGS 

Three  different  loading  configurations  were  chosen  for  Case  2, 
uniform  tension,  three-point  bending  with  S/W  of  4  and  8  (Figure  20). 

Each  of  the  different  Case  2  loading  configurations  requires  another 
function  p(t)  as  the  noncracked  body  stress  at  the  crack  location  stresses. 


a°  UNIFORM  TENSION  °° 

_ i _ 

==E 

THREE-POINT  BENDING 

_ 1 _ 

L  _ 

‘w 

P, 

'2 

/2r  s  n p  > 

Figure  20.  Case  2  Loading  Configurations 

5.  UNIFORM  TENSION 
In  this  case. 


p(t)  =  const.  =  o 

o 


(37) 


where  a  is  the  remote  stress.  Since  the  standard  nondimensional 
o 

representation  for  Kj  for  uniform  tension  Is 

K 

TTaa  <38) 

o 

the  function  p(t)  should  be  factored  by  l//rraa0  to  get  Kj  in  the 
desired  form.  Therefore,  take 


Pit)  =  — 

/rra 


(39) 


(in  the  computer  program  l/Zira  is  already  factored  in  the  weight  function 
so  p(t)  is  substituted  as  1.). 
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The  results  of  the  nondimensional  stress  intensity  factor  for 
uniform  tension  at  equally  spaced  points  from  a/W=.l  to  a/w=.7  are 
presented  in  Table  4.  These  results  are  compared  to  Tada's  (Reference  12) 
solution.  The  results  agree  with  the  reference  solution  within  8  percent. 

Figure  21  compares  the  reference  solution  (Reference  12)  to  the  values 
calculated  at  the  same  a/W  points  where  experiments  were  performed  on 
the  four-point  bending  specimen.  One  can  see  that  the  last  point  is  the 
only  calculated  value  that  does  not  agree  favorably  with  the  reference 
data.  Since  the  computations  required  differentation  of  experimental 
data,  one  should  be  very  suspicious  about  the  results  obtained  at  the 
experimental  interval  end  points.  Close  to  the  end  points  computing 
derivatives  actually  requires  extrapolations  that  can  be  very  inaccurate, 
since  the  behavior  of  the  function  is  unknown  in  these  regions. 

Investigating  Equation  30  shows  clearly  that  the  derivatives  3K/3a 
and  3nQ/3a  have  a  relatively  larger  influence  on  the  value  of  the  weight 
function  at  longer  crack  lengths;  therefore,  the  computed  values  are 
more  likely  to  get  distorted  at  the  upper  limit  of  the  experimental  a/W 
interval  as  can  indeed  be  seen  in  Figure  21. 

6.  THREE-POINT  BENDING 

A  solution  for  the  stress  distribution  of  a  noncracked  strip  under 
three-point  bending  in  a  cross-section  through  the  center  is  given  by 
Timoshenko  (Reference  27): 


o 


x 


+  2ircb 


A  i\ 

10  c) 


(40) 
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TABLE  4 

Kj  (NONDIMENSIONAL)  FOR  UNIFORM  TENSION 


Crack  Length 

a/W 

KI  .... 

Relative 
Difference 
in  % 

3  X  C  U  X d  u  CQ 

o/tT a 

Using  the 
Weight  Function 

“  r  x.  uin 

o/ira 

Ref  12 

.1000 

1.2705 

1.1957 

6.26 

.1500 

1.2697 

1. 2682 

.12 

.  2000 

1.4041 

1.3667 

2.74 

.  2500 

1.5540 

1.4941 

4.01 

.3000 

1.6855 

1.6551 

1.83 

.3500 

1.8170 

1 

1.8565 

-2.13 

.4000 

1.9967 

2.1080 

00 

CM 

in 

1 

.4500 

2.2917 

2.4241 

-5.  46 

.  5000 

2.7719 

2.8266 

-1.93 

.  5500 

3.4752 

3.3486 

3.78 

.6000 

4.3636 

4.0432 

7.92 

.6500 

5.3158 

4.9993 

6.33 

.7000 

5.1832 

6.3755 

-3.02 

I 
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The  terms  of  this  equation  are  defined  in  the  following  Figure  22. 
ox  is  the  stress  in  x-direction  along  AD. 


Figure  22.  A  Uniform  Strip  Under  Three-Point  Bending 
According  to  Timoshenko,  this  expression  gives  the  stress  with  very  good 

p 

accuracy  except  for  point  D  where  an  error  of  .121  ^  exists  based  on  a 
more  accurate  solution.  This  term  will  lead  in  our  worst  case  to  an 
error  of  4  percent  in  the  load  distribution  at  the  crack  mouth  location, 
but  since,  at  this  point  the  weight  function  has  its  lowest  value,  the 
integral  of  the  product  of  h  and  p(t)  is  expected  to  have  much  better 
accuracy.  The  standard  nondimensional  solution  is  of  the  form 


a  /ifa 


(40a) 


where 


6M 

o  =  - 

4  be 

So  the  undimensional  form  of  Kj 


-  and 

> 

becomes 

Kj  4c2b 
/iTa  3P£ 


M  = 


EA 

2 


(41) 


40 
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The  factor  l//rra  is  already  included  in  our  computer  program  in  the 
weight  function  so  the  form  of  p(t)  that  is  required  in  order  to  get  Kj 
in  the  form  of  Equation  41  is 

o  (t) 4bc2 

=  3?I  (42) 


Substituting  Equation  40  to  Equation  42  yields 


where  S  =  2£ 
and  W  =  2c 


(43) 


Changing  to  crack  tip  coordinates  t  the  expression  for  y/c  becomes 

y/c  =  7w72r(fc+f  _a)  =  2  (S  -  (?) + 1  (4 


Substitute  Equation  44  into  Equation  43  and  get 


(45) 

Results  for  the  nondimensional  stress  Intensity  factor  considering 
three-point  bending  (S/W=4  and  S/W=8)  are  presented  in  Tables  5  and  6 
respectively.  These  results  are  compared  with  known  results  from 
References  12  and  28.  The  reference  solution  claims  to  be  accurate  for 
a/W  <.6  and  In  this  interval  the  calculated  experimental  values  are 
within  6.15  percent  accuracy  for  S/W=4  and  5.9  percent  for  S/W=8. 

Figures  23  and  24  show  the  calculated  Kj  for  three-point  bending  for 
the  a/W  points  where  actual  data  was  gathered  for  the  four-point  bending. 

As  for  the  uniform  tension  case,  the  last  point  Is  the  only  one  not  fitting 
the  expected  curve  because  of  the  inaccurate  extrapolation  procedure 
explained  previously. 
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TABLE  5 

Kj  (NONDIMENSIONAL)  FOR  THREE-POINT  BENDING,  S/W=4 


Crack  Length 

a/W 

(Nondiinensional) 

Calculated  Using 
the  Weight  Function 

Kj  (Nondimensicnal) 
Ref  28 

Relative 
Difference 
in  % 

•  1000 

1.0134 

.9849 

2.89 

.1500 

.9447 

.9731 

-2.91 

.2000 

.  9838 

.9803 

.  36 

.  2500 

1.0272 

1.0037 

2.34 

.3000 

1.0491 

1.0425 

.64 

.3500 

1.0635 

1.0979 

-3.13 

.4000 

1.1015 

1.1736 

-6.15 

.4500 

1.1994 

1.2754 

-5.96 

.5000 

1.3866 

1.4112 

-1.74 

.  5500 

1.6656 

1.5910 

4.69 
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TABLE  6 

Kj  ( NONDIMENS IONAL )  FOR  THREE-POINT  BENDING,  S/W=8 


Crack  Length 

a/W 

Kj  (Nondimensional) 

Calculated  Using 
the  Weight  Function 

Kr  (Nondirnensional) 

Ref  28 

Relative 
Difference 
in  % 

.1000 

1.0630 

1.0156 

4.66 

.1500 

.9911 

1.0080 

-1.68 

.2000 

1.0314 

1.0183 

1.28 

.2500 

1.0755 

1.0438 

3.04 

.  3000 

1.0964 

1.0839 

1.15 

.3500 

1. 1087 

1.1403 

r- 

CM 

l 

.4000 

1.1450 

1.2167 

-5.89 

.4500 

1.2431 

1.3189 

-5.75 

.  5000 

1.4328 

1.4549 

-1.52 

.5500 

1.7160 

1.6349 

4.96 

43 
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SECTION  VI 
CONCLUSIONS 

1.  This  work  shows  clearly  that  the  experimental  weight  function 
procedure  leads  to  excellent  results  for  the  stress  intensity  factor 
without  actually  testing  the  loading  configuration  in  question. 

2.  A  set  of  experiments  of  a  single  loading  configuration  can  be 
used  to  calculate  the  stress  intensity  factor  solution  for  any  other 
loading  configuration  applied  to  the  same  geometry.  The  only  additional 
information  required  is  the  stress  distribution  at  the  crack  location 
for  a  noncracked  body. 

3.  One  should  limit  the  results  obtained  for  a  different  loading 
only  to  the  same  a/W  region  studied  in  the  original  experiment.  Extra¬ 
polations,  or  even  calculations,  near  the  test  a/W  limits  may  lead  to 
inaccurate  results. 

4.  Because  data  scatter  can  lead  to  errors  when  differentiated,  care 
should  be  taken  to  collect  enough  data  to  yield  reasonable  interpolating 
functions. 

5.  Even  though  some  error  was  introduced  in  the  Kj  and  n  measure¬ 
ments,  the  numerical  procedure  using  the  derivatives  of  Kj  and  nQ  did  not 
distort  the  Case  2  results  beyond  the  basic  accuracy  of  the  Case  I  measure¬ 
ments. 

6.  The  laser  interferometric  technique  provides  an  efficient  procedure 
for  crack  opening  measurements  and  leads  to  accurate  results  for  this  type 
of  measurement. 

7.  The  Orange  (Reference  11)  conic  section  again  provided  an  accurate 
representation  for  the  crack  opening  profile  for  edge  cracked  strips.  This 
conic  section  should,  however,  be  compared  with  actual  profile  measurements 
(as  done  here)  prior  to  applying  it  to  other  geometries. 

8.  The  fact  that  the  actual  cracks  were  not  perfectly  straight  on  a 
micro-scale,  as  assumed  in  this  theory,  is  the  most  likely  source  for  Kj 
calibration  errors  for  Case  1  loading  (four-point  bending).  Perhaps  other 
materials  would  give  "straighter”  cracks  which  more  closely  match  the 
analytical  behavior. 
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APPENDIX  A 
THE  WEIGHT  FUNCTION 

The  idea  of  the  weight  function  has  been  introduced  by  Bueckner 
(Reference  2)  and  was  also  discussed  in  References  2  through  8. 

The  following  analysis  shows  that  If  the  complete  solution  of  the 
stress  intensity  factor  and  the  crack  opening  displacement  for  a  crack 
subjected  to  a  certain  loading  system  Is  known,  then  the  stress  Intensity 
factor  solution  for  the  geometry  under  another  loading  configuration  may 
be  obtained  directly  from  the  known  solution. 

Let  us  consider  a  cracked  body  with  loads  P^,...,P  as  described 
in  Figure  A-l . 


Figure  A-l.  Loaded  Cracked  Body 
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The  elastic  energy  release  rate,  often  referred  to  as  the  Griffith  energy 
rate,  can  be  defined  following  Irwin  (Reference  22)  as 


'P 


(A-l ) 


where  U  is  the  total  potential  energy  stored  in  the  cracked  body  and  a  is 
the  crack  length. 


g  has  units  of  force  and  is  often  referred  to  as  the  crack 
extension  force 


U 


l  N 

7  I 


N 


i=l  j=l 


C  .  .P  .P  . 
13  i  3 


(A -2) 


where  C^j  is  the  compliance  coefficient  defined  as  the  deflection  at  point 
i  due  to  a  unit  force  at  point  j. 


Q=  1U 
c y  3a 


P 


N  3C .  .  ,  N  5u, 

X  H  P  P  =  I  X  i 

.  .  >  3a  l  j  2._.  3a 

1=1  g=l  J  1=1 


1  N 
k  I 


P  . 
1 


( A-3) 


where 


N  N 

l.  =  l  u.3  =  l  C .  .p  . 
1  j=l  1  j=l  *3  3 


( A- 4 ) 


UiJ  is  the  displacement  at  i  due  to  load  Pj. 

Since  the  elastic  energy  release  rate  g  (based  on  strain  energy 
consideration)  and  the  stress  Intensity  factor  K  (based  on  crack  tip 
stress  considerations)  are  both  parameters  that  characterize  crack  growth, 
a  relation  can  be  established  between  them  as  shown  by  Hertzberg 
(Reference  15). 


<5  =  IT  (A'5) 


48 


AFWAL-TR-80-4001 


where  H  is  defined  as 


H 


E 


E 

1-v2 


for  plane  stress 
for  plane  strain 


(A-6) 


Since  K  is  the  linear  elastic  stress  Intensity  factor,  superposition 
can  be  used  to  write 


K 


N 


i-1 


K. 

1 


N 

I 

i=i 


k. 


1 


(a^ 


(A-7) 


k.  stress  intensity  factor  per  unit  load.  Substituting,  into  the  elastic 
energy  release  rate 


N 

l 

i=l  j=l 


(a)  kj  (a) 


( A-8) 


By  equating  Equation  A-8  to  Equation  A-3  we  conclude  that  they  should  be 
equal  term  by  term. 

Thus, 


(a)  k^  (a)  x  3Ci^  (a) 
H  =  2  3a 


(A-9) 


If  the  full  solution  for  K  is  known  only  for  one  load,  say  P  „,  then 

m 


ki(a) 


H 

2 


3C.  (a) 

lm 

3a 


(A-10) 
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recall,  that  we  defined  in  Equation  A-7  k 


m 


ui 

and  in  Equation  A-4  C^m  =  p- 


uj_ 

m 


Therefore, 


~  m 

k  ,a)  .  a  X 

i,a)  2  3a  K 
m 


(A-ll) 


and  since 


N 


K  =  y  k .  (a)  P  . 

•  “  i  l  i 

i=l 


(A- 12) 


Then 


u  N  3  u  -m 

v  —  _JL_  y  _ p 

K  2K  .  i1  3a  i 
m  i=l 


(A-13) 


For  an  arbitrary  distribution  of  ,  and  a  set  of  surface  tractions  T(s), 
the  expression  for  K  can  be  written  as 


K  _  _iL_  f  Of 

K  “  2K  3a 

m  J 

*'  C? 


(s,m)  T(s)ds 


( A- 14) 


If  one  defines  the  weight  function  as 


H _  3  um  ( s ,  a ) 

9  cl 


hm(s,a) 


2Kla) 

m 


(A-15) 


Then  we  can  readily  see  that  the  weight  function  depends  only  on  the 
geometry  of  the  cracked  body  (Including  the  crack  length)  and  the 
loading  configuration  denoted  by  the  superscript  m. 
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Substituting  the  term  for  the  weight  function  Into  Equation  A- 14, 
we  get 


K 


h  (s,a)T(s)ds 
m 


(A- 16) 


We  can  therefore  conclude  that  the  stress  Intensity  factor  for  the 
arbitrary  surface  traction  T(s)  can  be  obtained  from  the  known  weight 
function  we  computed  considering  the  loading  configuration  m.  In  a 
little  bit  more  general  derivation,  the  body  forces  can  be  included  and 
one  can  write 


K 


T-hdr 


f • h*  dA 


(A-1  7) 


where  h  is  the  weight  function  obtained  from  any  known  load  configuration, 

T  are  the  surface  traction  for  which  K  is  to  be  determined, 

f  are  the  body  forces  for  which  K  is  to  be  determined, 

[’  is  a  chosen  path  around  our  specimen  as  to  include  all  the 
surface  tractions  T,  and 

A  is  the  region  where  the  body  forces  act. 


SI 
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APPENDIX  B 

JUSTIFICATION  OF  ASSUMPTIONS 

1.  USE  OF  LINEAR  ELASTIC  FRACTURE  MECHANICS 

Linear  elastic  fracture  mechanics  (LEFM)  requires  that  the  amount 
of  plasticity  near  the  crack  tip  will  be  relatively  small.  As  a  "rule 
of  thumb"  LEFM  can  be  considered  to  be  valid  if 

rp  <  0.1a 

where  rp  is  the  radius  of  the  plastic  zone  around  the  crack  tip  and  a  is 
the  crack  length.  Using  Irwin's  circular  model  for  the  plastic  zone 
incorporating  the  upper  bound  dimension  as  obtained  by  a  plane  stress 
solution  (Figure  B-l ) , 


Figure  B-l.  Irwin  Circular  Plastic  Zone  Model 


where  o  . $  is  the  yield  strength  of  the  material. 
Using  7075  T651 

°y's  =  80.8  ksi 


(B-l) 
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The  maximum  Kj  used  in  our  tests  was  lower  than  15  ksi/in. 


Therefore, 

the  maximum  radius  of  plastic  zone  that  has  been  developed  was 

2 


=  _L  =  0.0055  in 

2tt  \80. 8/ 


The  smallest  crack  length  used  was 

a  =  0.0815" 


Thus , 


max 


0.0674  <  0.1 


The  conclusion  is  that  LEFM  was  valid  even  in  the  most  extreme  conditions 
In  most  cases  the  crack  length  was  significantly  longer  and  therefore 
rp/a  smaller. 


2.  PLANE  STRESS/STRAIN  CONSIDERATIONS 

As  a  rule  of  thumb  (References  15,  16),  if 

/Kto  '2 
b 


1 2'5fe) 


(B-2) 


where  b  is  the  material 
the  yield  strength. 


thickness,  Kj^ 


the  fracture  toughness,  and  o 


y's 


The  specimen  can  be  considered  thick  enough  to  be  under  a  plane 
strain  stress  field. 


Kjc  for  7075  T651  plate  is  21.4  ks i >/Tn . 


2.5 


0. 176" 


Since  .24"-  .25"  thick  specimens  were  used,  they  can  be  viewed  to 
be  under  plane  strain  conditions. 
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APPENDIX  C 

LASER  INTERFEROMETER  DATA  REDUCTION 

A  typical  strip  chart  recorder  trace  of  fringe  motion  crack  mouth 
opening  and  load  is  shown  in  Figure  C-l.  From  this  chart  the  fringe 
order  versus  load  data  can  be  obtained  for  the  left  and  right  fringe 
pattern  as  presented  in  Table  C-l.  The  data  can  be  plotted  as  shown 
in  Figure  C-2. 

From  Figure  C-2  slopes  of  Sm/P  of  the  left  and  right  fringe  pattern 
can  be  obtained  by  considering  only  the  linear  part  and  ignoring  the 
initial  nonlinear  part  of  the  plot  that  is  associated  with  fatigue  crack 
closure  effects.  Those  slopes  are  averaged  to  eliminate  free  body 
motion. 

The  experiment  is  repeated  and  results  are  averaged.  The  value  of 
the  slope  is  used  in  Equation  33  to  get  displacement  per  unit  load  data. 
Through  Equation  34  stress  intensity  factors  can  be  calculated. 


Figure  C-l.  Typical  Stripchart  Recorder  Trace  of  Fringe  Motion,  Load 
and  Crack  Mouth  Opening 
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APPENDIX  D 

COMPUTER  PROGRAM  DESCRIPTION 

The  computer  program  has  been  written  to  perform  the  numerical 
computation  involved  in  the  process  of  data  interpolation,  weight  function 
calculations  and  finally  evaluation  of  the  integral  in  Equation  30  to 
determine  for  Case  2. 

In  order  to  allow  maximum  flexibility  in  the  usage  of  the  program, 
it  was  constructed  of  a  main  program  that  performs  the  integration  and 
nine  subprograms  that  perform  all  the  other  necessary  calculations.  In 
addition,  two  library  subprograms  are  used,  one  from  the  IMSL  library 
and  the  other  from  the  CC  6600  library.  This  structure  of  the  program 
allows  using  it  for  different  cases  and  configurations  and  permits  the 
user  to  enter  his  own  interpolation  schemes,  comparison  functions  and 
the  noncracked  body  stress  distribution  p(t)  that  makes  it  possible  to 
compute  any  desired  Case  2  stress  intensity  factor. 

The  program  is  written  in  FORTRAN  and  was  executed  on  the  CDC  system 
at  Wright-Patterson  Air  Force  Base.  Execution  time  for  the  three-point 
bending  and  uniform  tension  Case  2  stress  intensity  factor  varied  from 
two  to  four  seconds. 

1.  THE  MAIN  PROGRAM 

The  main  program  reads  in  the  input  data,  calls  in  the  different 
subroutines  to  perform  necessary  computation  and  performs  the  numeric 
integration  by  using  the  IMSL  library  subprogram  DCADRE.  An  iterative 
scheme  "pushes"  down  the  lower  integration  limit  until  the  desired  accuracy 
for  the  end  result  is  met.  At  the  end,  it  prints  out  computed  Kj  for 
Case  2  compared  to  user-supplied  comparison  data. 

2.  FUNCTION  SUBPROGRAM  F(x) 

This  function  computes  the  value  of  the  integrant  (the  weight 
function  multiplied  by  P(t))  at  each  distance  t  from  the  crack  tip. 
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3.  FUNCTION  SUBPROGRAM  P(x) 

This  function  computes  the  value  of  the  stress  at  the  crack  location 
P(t)  for  a  noncracked  body  at  any  distance  t  from  the  corresponding  crack 
tip. 


The  function  should  be  supplied  by  the  user  and  is  the  only  subprogram 
that  has  to  be  altered  when  computing  Kj  for  different  Case  2  configurations. 
Care  should  be  taken  as  related  to  the  form  in  which  P(t)  is  supplied. 

This  form  will  determine  the  form  of  the  calculated  Kj.  A  factor  of 
l//ii¥  is  already  included  In  the  weight  function  calculation.  So  if 
P(t)  is  supplied  directly  in  terms  of  stress,  the  resulting  Kj  values  will 
be  in  the  dimensional  form  of  Kj/Zira. 

4.  LIBRARY  FUNCTION  PLSCF 

A  CC  6600  Library  subroutine,  this  function  computes  least  square 
polynomial  fits  for  a  given  set  of  data.  Polynomials  up  to  6  degrees  or 
Chebyshev  polynomials  of  any  degree  can  be  specified  by  the  user.  For 
further  information  see  Reference  25. 

5.  LIBRARY  FUNCTION  SUBPROGRAM  DCADRE 

An  International  Mathematics  and  Statistics  IMSL  Library  subprogram, 
this  subprogram  performs  numeric  integration  of  a  function  using  the 
Romberg  extrapolation  scheme.  The  function  DCADRE  is  capable  of  integrating 
functions  with  "jump"  discontinuities  and  certain  types  of  singularity. 

This  feature  is  most  important  in  our  case  because  of  the  crack  tip  singu¬ 
larity  of  the  weight  function.  Even  if  the  desired  accuracy  is  not 
found,  DCADRE  returns  the  best  available  estimate.  See  the  IMSL  Manual 
(Reference  26)  for  further  information. 

6.  SUBROUTINE  KPOL 

This  subroutine  computes  a  fourth  order  least  square  polynomial  fit 
for  the  nondimensional  Kj  data  using  PLSCF  subroutine  from  the  CC  6600 
library.  It  prints  out  the  polynomial  coefficients  C(5)...,C(1)  and 
returns  these  values  to  the  main  program.  It  also  prints  out  the  error 
flag  IER  from  PLSCF. 
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7.  SUBROUTINE  EPOL 

This  subroutine  computes  a  fourth  order  least  square  polynomial  fit 
for  the  nondimensional  n0  data  using  PLSCF  subroutine  from  the  CC  6600 
library.  It  prints  out  the  polynomial  coefficients  D(5) , . . . ,D(1 )  and 
returns  them  to  the  main  program.  It  also  prints  out  the  error  flag 
IER  from  PLSCF. 


8.  SUBROUTINE  EVEK 


Evaluates  the  nondimensional  Kj  from  the  fourth  order  least  square 
polynomial  and  the  derivative  of  Kj  with  respect  to  crack  length  3Kj/3a 
at  the  points  where  Kj  for  Case  2  should  be  evaluated.  One  should  realize 
that  the  differentiation  considers  the  fact  that  the  nondimensional  Kj 
usually  factors  out  1/vci —  that  should  be  included  in  the  differentiated 
term.  The  derivative  3Kj/3a  requires  the  expression 


5KI  _  3  (o/TTa  Y) 

3  ci  3^ 


where  Y  is  the  nondimensional  Kj  value  and  a  is  the  maximum  tensile 
stress  in  the  specimen.  Kj  and  3Kj/3a  are  returned  to  the  main  program 
in  a  dimensional  form  factored  only  by  1/a  as  to  match  the  n  data. 

9.  SUBROUTINE  EVETTA 

Evaluates  the  nondimensional  nQ  from  the  least  square  fourth  order 
polynomial  and  the  derivative  of  nQ  with  respect  to  crack  length  3nQ/3a, 
at  the  points  where  Kj  for  Case  2  should  be  evaluated.  One  should  realize 
that  the  differentiation  considers  the  fact  that  the  nondimensional  n 

o 

is  factored  by  a  term  of  "a"  that  should  be  included  in  the  differentiated 
term  as  applied  in  the  above  paragraph. 

nQ  and  3g  /3a  are  returned  in  a  dimensional  form  factored  by  1/a 
to  match  the  Kj  data. 
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10.  SUBROUTINE  COPK 

Compares  the  nondimensional  Kj  input  data  with  the  interpolated  data 
by  KPOL  at  the  test  points.  A  comparison  table  that  includes  a  relative 
error  in  percents  is  printed  out  for  the  user's  evaluation  and  judgement. 

11.  SUBROUTINE  COPETA 

Compares  the  nondimensional  nQ  input  data  with  interpolated  data 
by  EPOL  at  the  test  points.  A  comparison  table  that  includes  a  relative 
error  in  percents  is  printed  out  for  the  user's  evaluation  and  judgement. 

12.  SUBROUTINE  KANA 

Evaluates  Kj  for  the  Case  2  loading  configuration  from  a  user- 
supplied  comparison  function.  The  evaluation  is  done  at  the  same  points 
where  Kj  Case  2  is  evaluated  from  the  experimental  data.  This  subroutine 
can  be  eliminated  if  the  above  comparison  data  is  available  explicitly. 

In  this  case,  the  comparison  values  should  be  assigned  directly  to  AK 
array.  If  no  comparison  values  are  available,  the  final  printout  table 
should  be  altered  to  avoid  undefined  terms. 

13.  LIST  OF  SYMBOLS  USED  IN  THE  PROGRAM  (ALPHABETICAL  ORDER) 

A  Constant  used  in  subroutine  KANA 

AERR  Input,  absolute  error  in  the  integral  routine  DCADRE 

AK(99)  Kt  /a/i ta  comparison  value  of  the  KT  Case  2 

lQ)  1 

nondimensional  evaluated  from  the  user-supplied 
comparison  function  at  points  defined  by  S 

AKS ( 99 )  Kt  /a/na  calculated  values  of  K,  Case  2 

'©  1 

(nondimensional)  at  points  defined  by  S 

B  Interim  term  for  weight  function  calculations 

C(5)  Interpolating  polynomial  coefficients  for  Kj  data 

D(5)  Interpolating  polynomial  coefficients  for  nQ  data 

DB  Interim  term  for  weight  function  calculation 

DC  3R/3a  interim  term  for  weight  function  calculation 
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| 

1 

l 

1 

? 

DY(99) 

Interim  derivative  term  for  DYA 

DYA ( 99 ) 

1 /a  (3Kj3a)  derivative  w*r*t  the  crack  length  of  Kj  j 

Case  1  at  points  defined  by  S  ( 

D2(99) 

Interim  derivative  term  for  DZA 

DZA ( 99 ) 

l/o  (3nQ/3a)  derivative  rrt  the  crack  length  of  nQ  , 

at  points  defined  by  S 

| 

EDS 

Input,  desired  absolute  accuracy  in  Kj  computation  | 

ES ( 99 ) 

(AKS-AK)  100/AK  relative  error  in  the  computed  Kj  1 

EST 

Current  estimate  of  value  of  the  integral  by  DCADRE 

E Y (99 ) 

(YC-Y)  1 00/ Y  relative  error  in  percents  of  interpolated  f 

Kj  data  at  points  defined  by  S 

EZ ( 99 ) 

i 

(ZC-Z)  100/Z  relative  error  in  percents  of  interpolated 
nQ  data  at  points  defined  by  S 

F 

h,  the  weight  function  (subfunction) 

G 

Upper  limit  of  the  integral  [ 

H 

l 

Input,  H  material  property  j 

IER 

Error  flag  for  library  subprograms  PLSCF  and  DCADRE  ! 

M 

Input,  number  of  experimental  data  points 

N 

Input,  number  of  points  where  KT  should  be  evaluated  for 

Case  2  1 

P 

Stress  distribution  for  the  noncracked  body  at  the  crack 
location  for  Case  2  loading  configuration  (subfunction) 

PI 

IT 

Q 

Lower  limit  of  integral  1 

R 

Interim  term  for  weight  function  calculation  < 

RELERR 

Input,  relative  error  in  the  Integration  routine  DCADRE  ! 

S  (99 ) 

Input,  a/w,  nondimensional  crack  length  at  point  where 
experiments  were  performed 

T 

Interim  term  for  weight  function  calculations 

V 

Previous  interim  value  for  EST 
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Card 

Variable 

No,  of 
Variables 

Real/  Format 

Number 

per  card 

Integer 

Real  Unformatted 
Real  Unformatted 
Integer  Unformatted 
Real  Unformatted 
Integer  Unformatted 
Real  Unformatted 
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t  i 

l  , 

l" 

VA 

(EST-V),  difference  between  last  and  current  Integration  j4 

estimate  ^ 

. 

W  ( 99 ) 

Weight  corresponding  to  each  data  point  used  in  PLSCF  ' 
(currently  defined  as  -1  to  give  equal  weight  to  each  [i 
point)  rj 

W0RK(25) 

Work  array  for  PLSCF  subroutine 

; 

X  ( 99 ) 

Input,  a/w,  nondimensional  crack  length  for  points  ! 

where  Kj  Case  2  should  be  evaluated 

\ 

Y  (99) 

Input,  Kj /  o/ira,  nondimensional  Kj,  experimental  data 

t 

YA(99) 

KI/a,  interpolated  Kj  for  Case  1  at  points  defined  by  x  . 

YC ( 99 ) 

Kj/  o/ira,  interpolated  nondimensional  Kj  for  Case  1  at  , 

points  defined  by  S 

YK(99) 

Kj/  o/i ref,  interpolated  nondimensional  Kj  for  Case  1  at 

points  defined  by  x  K 

i 

Z(99) 

Input,  nQ*H/(2aa),  nondimensional  crack  mouth  opening  j 

ZA(99) 

i 

n Jo,  interpolated  nQ  at  the  points  defined  by  x 

ZC(99) 

nQH/(2oa),  interpolated  nondimensional  at  points 

defined  by  S 

ZK( 99 ) 

! 

n0H/(2oa),  interpolated  nondimensional  at  points 

defined  by  x 

14.  SEQUENCE  OF 

, - 

INPUT  DATA 

-  .  : =  :  r  =z=i 
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15.  SEQUENCE  OF  OUTPUT  DATA 

1.  Interpolation  polynomial  coefficients  for  Kj 

2.  Interpolation  polynomial  coefficients  for  nQ 

3.  Comparison  table  of  actual  to  Interpolated  Kj 

4.  Comparison  table  of  actual  to  Interpolated  nQ 

5.  The  following  printout  will  appear  only  if  singularities  have 

been  met  by  subroutine  DCADRE.  Each  time  a  singularity  was  met,  a  warning 
statement  is  printed  out  followed  by  the  values  of  the  crack  length 

(X(I)),  the  distance  from  crack  tip  (Q),  and  the  value  of  the  error  flag 

(IER)  for  which  the  singularity  occurred.  Values  of  IER  equal  to  65  or 

66  indicate  that  the  singularity  was  successfully  handled.  If  any  other 
value  appears  consult  the  IMSL  Manual  (Reference  26). 

6.  Comparison  table  of  Kj  for  Case  2  calculated  values  and  supplied 
comparison  values. 

16.  LISTING  AND  RESULT  SAMPLE 

The  following  pages  include  a  complete  listing  of  the  program  and  a 
sample  printout  of  results  for  uniform  tension  as  Case  2  loading  config¬ 
uration. 
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ta/ta  opt*i 


FTN  *»  .  7  ♦*»  7  b 


I9/2A/79  1M7.I7  PAGE  | 


1 


9 


II 


19 


2 1 


29 


SI 


39 


M 


b9 
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PROGRAM  TH  5  (INPUWIO  .OUTPUT) 

C 

. . ••***•*•*  OCCCNBEK  1979 

C  DATA  REDUCTION  FOR  THSIS  -OAM  0AT-TIKVA 

C  (SEE  THS1S  APPENDIX  ( AFI T/GAE/AA/790-2  »  'OR  OCTAUEO  DESCRIPTION) 

C 

C 

01  HEN  SION  X  (99 ) , Y (99) » v (99) ( 5(991 , Y< t 991 , f Al 99) .VC (99)  ,0YA(99) ,DY( 

*99) .E  Y<99) ,7* (99) ,?A(9«) ,7C(99) ,)7AI99l ,07(99) ,AA(99> ,AAS(99) , ES(9 
#9) ,C<5) ,0(5), WORK (791 , M ( 99) , E 7 (93  I 
COHHON/OAN/G#T.»,R,Oe,nC 
EXTERNAL  f 
PI  ><»,  •ATAN(1«0) 

C  FIRST  CARO 

C  READ  IN  H 

READ*  i  H 

C  SECCONO  CARO 

C  R£  AO  IN  THE  FOLLOING  TERMS  ON  TH;  SAME  CARO! UNFORMATTED) 

C  AERR-ABSOLUT  ERROR  IN  THE  INTEGRAL  RO U T I NE (OCA ORE > 

C  REllPR-RE LA(IVE  T RROR  IN  THE  INTEGRA!  ION  ROU II NC ( OCAORE ) 

C  CPS-OESIREO  A9S0LUIE  ACCURACY  IN  ANSWER 

RE  AO*  *CRR,REIERR,EPS 

C  NEXT  C*  O-NUMBER  OF  EXPERIMENTAL  TATA  POINTS 
READ*  i  M 

C  THE  NEXT  N  CAROS  CONTATNE  THE  EXP:  RIME  NT  AL  DATA  THREE  VALUES  ON  EACH 

C  CARO-A/N,K(NCN  OIHEnSIONAU  ,ETA(N5N  01  MENSIONAD-  (UNFORMATTED) 

DO  2  IMiM 

2  REAO*  ,S(I)»T(I)  ,7(1) 

C  NEXT  C A RO-N-NUiBEA  OF  POINTS  WHERE  K  SHOULD  BE  EVALUATED  FOR  CASE  2 
RE  AO*  ,N 

C  NEXT  N  CAROS-VALUTS  OF  A/M  AT  mhi:h  A  SHOULO  BE  EVALUATED  (ONE  VALUE  PER 

C  CARO-UNFORHATTEO) 

DO  17  IM.N 
17  REAO*  i  x  ( I  ) 

PRINT  1 

i  FORMAT  IlHltllX, "RESULTS") 

PRINT*,-  ■•«•■•••* 

PRINT*,**  ** 

PRINT*,**  “ 

PRINT*,**  ** 

PRINT*,**  - 

C  COMPUTE  A  LEAST  SQUIRE  FORTH  OROER  POLYNOMIAL  FOR  A  DATA  A  HO  PRINT  OUT 
C  POLYNOMIAL  COEf  F l Tl f  NTS  r(S),,,,CII)  ANO  ERROR  FLAG  OF  SUBROUTINE  PlSCF 
CALL  XPOL <S,Y,M,C,MORK,M) 

C  COMPUTE  A  LEAST  SQUARE  FORTH  OROER  POLYNOMIAL  FOR  ETA  DATA  ANO  PRINT  OUT 
C  POLYNOMIAL  COEFFITTENTS  0(5). ...DIM  And  ERROR  FLAG  OF  SUBROUTINE  PLSCF 
CALL  I POL ( S , 7  *  N , D , MOR  * , M I 

C  EVALUATE  *  A NO  JfS  ICRIVATIVES  FT OH  THE  LEAST  SQUARE  POLYNOMIAL 
CALL  fVEA (X,C,PI,N, Y  A  ,  P  Y  A  ,  Y  A  ,  0  Y ) 

C  EVALUATE  ETA  ANO  ITS  QFPIVATIVES  F EON  THE  LEAST  SQUARE  POLYNOMIAL 
CALL  EVEETA(X,0,h,n, *«, 02 A, 7 A ,0*1 

C  COMPARE  THE  ACTUAL  EXPERIMENTAL  A  PAT  A  Ml  TH  THE  LEAST  SQUARE  FIT 
CALL  COP A (S,V,C,H,vC,FV) 

C  COMPARE  THE  ACTUAL  EX»EPlMENlAl  ETI  OATA  M l T  H  THE  LEAST  SQUARE  FIT 
CALL  COPETA($,  7 #0,M, *c,e/> 

C  COMPUTE  COMPARISON  VALUF*  FOR  A  IIS!  2  (USER  SUPPLIED  ROUTINE) 

CALL  AANA(X,N,AK,PI) 
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mograh  tm* 


74/7%  OPT*l 


FT  W  4.7*476 


09/20/79  11.07.02 


PRINT*, “  - 
PRINT*,4*  “ 

PRINT*, M  " 

PRINT*,**  - 

COMPUTE  K  FOR  CASE  2  AND  PRINT  OUT  A  COMPARISON  TABLE  BETWEEN  THE  CALCULATED 
VALUES  ANO  THE  SUPPLIED  COMPARISON  VALUES 

PRINT*, “THE  FOLLOING  PRINT  OUT  IS  FROM  THE  NUMERICAL  INTEGRATION** 

PRINT*, “SUBROUTINE  (OCADRE*  AND  UDICATS  THAT  SINGULARITIES  WHERE** 

PR  I  NT*,  “HANOI  E  O'* 

PRINT*, “  - 

PRINT*, “  ICR* ERROR  FLAG, IF  65  3R  6b  -SINGULARITIES  WHERE  SUCCESS 

•  FULLV  HANOL£D,ANV  ATHFP  VALUE  CONSULT  IMSL  HANUAL4* 

PRINT*, “  MXI-IHC  CRACK  LENGTH  FOR  WHICH  THE  ERROR  FLAG  HER)  IH 

•  PL  IfS“ 

PRINT*, “  O-OISTANCE  TROH  CRACK  TIP  RELATIVE  TO  WIOTM  (T/WI  FOR 

•  WHICH  THE  SINGULARITY  WAS  OfeTACTTO** 

PRINT*,"  “ 

PR  INI  *,“  “ 

PRINT*, “  “ 

00  6  1*1, N 

PRINT •  - 

PRINT  *,*•** 

B*  0 . • TtA( II)**27(PI*H**ZI 
08«lfc ,*YA  CII*OYA (It f (PI*M**2» 

R*  I7A  (I>/X(II  >**2-0.*  <  YA  <IM**2/(Pl*X  tI>fH«*2» 

0C*-2.MZA(IM**2/<*Ul)**l-l6.*rA<I)  /(PI*xm*H**2I*OfACI>*0.*UA 

•  (i  n**2/ 1 pi* in* xii >  i**?>  *2.*zaii)  *oz a  a)  /m in*** 

G* x  II  ) 

T«H/M  (PI*xm  )••  ,f)*VA(in 
V»  0. 

0*  .000001 
0*0*.  1 

1ST*  OCADRF (F,0,G, AERR, REIT RR, ERROR, IERI 

IF < If  R « EQ • 01  GO  TO  21 

PR  1NT*,“XI“, I, **!■**, XTtl 

PRINT*  Q 

PRINT  *,*‘IER*-,IER 

CONTI  HUE 

VA*EST-V 

IF (VA  .IE. EPS)  GO  TO  16 

V*£ST 

GO  TO  15 

CONTINUE 

AKSU  »  *EST 

CSIII  «  (AKSUI-AKdl  I*  1O0.FAKTI) 

CONTI  NUE 
PRINT*, “  “ 

PRINT  *,“  “ 

PRINT*, “  “ 

PRINT*,"  “ 

PRINT  *,"  “ 

PRINT*, “  “ 

PRINT*,"  “ 

PRINT*, *’  RESULTS  FOR  K  UNIFORM  TENSION  USING  THE  WEIGHT  FUNCTIO 

•  N  HE  T  HOO  I  K  -NON  01  Nr  NS  T  ONI  L** 

PRINT  •,** . . . . . . . 
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*ROCRAM  rws  74/74  opt.i 


FTH  4.7*479  99/20/79  19.07.02  PACE  ] 


119 


129 


929 


1  SI 


PRINT*,"  - 
PRINT*,-  - 

PRINT*,”  A/N  -  CRICK  LENGTH  RE.1TIVE  TO  WOTM  - 

PRINT*,-  K-CXP  -  VALUES  CALCULATED  THROUGH  THE  NEIGH  FUNCTION- 

PRINT*,-  K-CAL  -  SUPPLIED  COMP*  RISON  VAIUES" 

PR1NT*,-  C  -  RELATIVE  ERROl  IN  X“ 

PR INI • " 

PRINT*,-  - 

PRINT*,-  A/M  K-E7P  R-C1L  f- 

PRINT*,-  “ 

00  13  1*1 tN 

PRINT  10,1(11  ,AKSIII  ,AKIII,ES(n 
19  FORMAT  UX.1F9. 4, %X,  IF# . 4 , EX , 1F0 • t , 2X , IFA  •  2 > 

1)  CONTINUE 

STOP-ENO  OP  PROGRAM- 
ENO 


FJNCTION  F 


74/74  0PT*1 


FTN  4.7*479  09/29/79  19.07.02  PAGE  1 


l  function  f no 

EXTERNAL  P 

COMHO N/OAN/6,r  ,9,R,0n,nC 

F*  (T/2.1*  <9*09*  X*0C*X»"  2*2. *R*X)M (9* X*R*x** 21 •* .5 > *P ( XI 
9  RETURN 

ENO 


FJNCTION  P 


74/74  OPT  *1 


FTN  4.7*479  09/20/79  10.07.02  PAGE 


l  FUNCTION  P<X) 

CO HMON/OAN/C, 1,9,9,09, DC 
P*l. 

RETURN 

9  ENO 
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subroutine 

i 

9 

II 

19 

79 

subroutine 

l 

9 

II 

19 

71 

SUBROUTINE 

1 

9 


«POt  7A/79  OPf»l 


PTN  %.7»*7fr 


09/70 /T9  11.0 /.I?  PAGE  1 


SUBROUTINE  KPOL (X,V ,N,C,MO*K.NI 

DIMENSION  Kill,  TUI  ,C(*>  , NORM  II,  Nil  I 

NO  EC*  4 

NP«H 

Will* 

iN*a 

CALL  PlSCFIK,V,H,NP,NDEG,NMAX,C,tN,XO,)(0,IIOR»(.IERI 
PRINT*,**  " 

PRINT*  ,-  " 

PRINT*,-  INTERPOLATION  POLINONlAL  FOR  K 

PRINT*,-  - 

PRINT*,-  - 

PRINT*  ,-C  (9»--,CI9) 

PRINT*, -C<A)»-,C(AI 
PRINT  ♦  ,-C  <31  *-.C(  3) 

PRINT*, -C(2»*",C<Z» 

PRINT*, -C  (II*-, Cli) 

PRINT*,-  * 

PRINT*, "IER-  ERROR  FLAG-SEE  ASO  lOHPUTER  SUBPROGRAM  LIBRARY  GUIDE* 

PRINT*, -ICR*-, IER 

RETURN 

ENO 

CPOL  7A/7A  OPT  *1  FTM  n.7 *476  19/78/79  18.0T.07  PACE  1 


SUBROUTINE  CPQLU  ,  7  »M  ,  0 ,  WORK ,  M> 

01  HENS  ION  XI II, ▼( II ,0(SI , WORK (II, Mill 

NO  EG*  A 

NP*H 

In*  o 

KllU-l, 

CALL  PLSCP(X,7,M,NP,NOEC,NNAX,0,(N,xn,XO,MORP,IERl 
PRINT*,-  “ 

PRINT*,-  - 

PRINT*,-  INTERPOLATION  POLINOHUL  FOR  ETA  - 

PRINT*,"  - 

PRINT*, **0(5)  *-,0(51 

PRINT*  ,"0 (A)*-, 0(A) 

PRINT*, -0(31  *-,0(  V 
PRINT*, -0(71*-, 0(7) 

PRINT*, -0(1)*-, 0(11 
PRINT*,-  - 

PRINT*, -IER-  ERROR  PL AC -SEE  ASO  JOMPUTER  SUBPROGRAM  LIBRARY  GUIDE- 

PRINT*  , -I  ER*-,  IER 

RETURN 

ENO 


EVER  7A/7A  0»*T«1 


ETN  A.MAT6 


09/78/79  18.07.17  PACE  1 


SUBROUTINE  EVEK (X ,C ,PI , N, TK , Or A , f «  ,  Of  I 
01  MEN  SION  Kill  ,C(9I»VR(1I,0VA(1I,VA(1),0YI1I 
00  1  1*1, N 

YK(II*C(5l*X(I)*,A*C(A)*XUI**J*:(T)*X(U*»7*C(?)*X(II*Cltl 
YA  II)  «YK(H*  (X(  I)  *PD  **  .5 

or  12)  »A*C<5)*XI  I)  **3*3'C(A)*XII)**7*?.*Cm*X(I)  *C  17) 

1  OVA  (I  |*CPI**.9l*(Yi(Hl/Cf  •**(1>I*0YCI1  l*XIII**l,9 

RETURN 
END 
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SUBTOUT XNE 

1 

s 

II 

If 

I« 

2f 

SUBROUTINE 

1 

f 

II 

if 

II 

If 


r 


\ 

i 

r 

t 

i. 

i' 


COPK  74/74  OPT-1 


FTN  4.7 *476 


09/11/79  11.07.02  PAGE  1 

I 


II 


S 

A 


SUBROUTINE  COPK(X,Y,C,M,YC,EY) 

01  HE  NS  ION  X<1>  ,Y<  1)  ,C  f  I  ,  VCU>  ,EM1> 

00  10  1*1 »H 

yc  ui  *c<* >**iz>  *#v*c  uipxcn”3»:i3>*x<i»*»*»cm*xui  ♦cm 


CY  III  ■  CYC  <I»-YUI  >*1|I./Y(I> 

PR 1NT  *  ,“  " 

PRINT*,**  “ 

PRINT*,**  “ 

PRINT*,**  ** 

PRINT*  ,**COHPERXSON  OF  ACTUAL  OATA  TO  INTERPOLATED  OATA* 
PRINT*,**  ** 

A/M  -  CRACK  LENGTH  Rt.ATIVC  TO  HIOTH  - 
K-EXP  ACTUAL  OATA" 

INTERPOLATED  DATA* 

RE  LA  T I  VC  ERROT  IN  X** 


K-CAl 

E 


PRINT  *,“ 

PR  1  NT  •  ,** 

PRINT  • 

PRINT  *,** 

PRINT*,-  " 

PRINT*,-  “ 

PRINT*,-  A/M  K-EXP  K-CAL  E  " 

PRINT*,-  - 
DO  A  1*1, H 

PRINT  3,X<I»,f <X> ,YCCI> ,EYCI> 

FORMAT  <lX,ir6.t ,AX,iri.%f2X»lF0.t,?X,  1FA.IT 

CONTINUE 

RETURN 

END 


COPEfA  74/74  Ot»T*l 


FTN  A. 7*476  19/26/79  16.07.02  PAGE  1 


19 


f 

7 


SUBROUTINE  C0PETA(X,7,0,M,7G ,€?> 

oi  hen  si  on  mii  ,  m»  ,cMr>  ,  zcu>  ,£r  m 

DO  19  1*1,  N 

?CtI)*0(5»*XCII**A*0(A)*X<n**3*)(3>*xa»**l*0(2)fxa»*0(l» 
E?U)  *  UC  (11-7  Tin*  ltO./Z<I> 

PRINT*,-  - 

print*,-  ** 

PRINT*,-  - 
PRINT* 

PRINT* 

PRINT* 

PR  INT  * 

PRINT  * 

PRINT* 

PRINT* 

PRINT* 

PRINT  * 

PRINT* 

PRINT  * 

PRINT* 

00  7  I»1,H 

print  ,rc m*E7<i> 

FORMAT  (1X,1FA.A,AX, lFt.A,2X,tF6.k,2X,lF6.2) 

CONTINUE 
RETURN 
END 


'  ,-COHPERISON  OF  ACTUAL  OAT1  TO  INTERPOLATED  DAT A* 

A/M  -  CRACK  LENGTH  RE.  ATI  VE  TO  WIDTH  - 
ETA-EXP  ACTUAL  OATA- 
,-  ETA-CAL  INTERPOLATED  JATA" 

(  RELATIVE  CRROE  IN  X* 


A/M  CTA-EXP  ETA-CIL 
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SUBROUTINE  EVEETA  TH/7k  0*»T »l 


FTN  4.7*479 


H/ll/M  tl.ir.12 


SUBROUTINE  EVEETA  (X  , 0  ,  H ,  N ,  ZK  , Oik,  ZA ,  0?) 

OZNENSZON  X(U,  0($>  ,  Z  X  (1  >  ,02  A  U  > ,  7  k  1 1 )  ,D?U> 

DO  5  I  *1 1  N 

ZK  <1)  •0m*KU)**4*0(4>*XUI**JO<IMX(I)**2*0(f>*X<IM0ti) 
ZA  (2WKI  X)*2.*XfI)/H 

OZ  (XI  •4*0(9)  *X(n  •*T>3*Dt4)*X<II**2*f  •*Ot$)*X(I>  ♦0I|> 

oz  a  t  j )  •  (t./H»*(om)*vm»ZK(oi 

RETURN 

(HO 


SUBROUTINE  KAMA 


FTN  4.7*474 


SUBROUTINE  KANA < X , N , A K ,PI ) 

OZNENSZON  XU)iAK(l) 

00  20  I-l.N 
A-X<J)*PZ/2. 

AK<Z)*(<2.*7AN(A>/(»Z*V(I>>>  **.f >*l  f • 792*2, ***X II) »• 37* tl.»SJNf A J 
•  )**S. }/CO$(AM 
CONTINUE 
RETURN 
END 
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RESULTS 


INTERPOLATION  POLINOH1AL  FOR  K 


C  (5)-7|. 75542913721 
C  (4)  — 92.45577  346502 
C(3)-42.9545S9M156 
C!21—  6.926729064317 
C  (I) *1.256007 10 54 34 

IER-  ERROR  FLAG-SEE  ASO  COMPUTER  SUBPROGRAM  LIBRARY  GUIDE 
I  CRM 


Interpolation  polinomiai  for  eta 

0(51-210. 3765959491 
0 16 )  —282*5727  32332 
0(31-150. 7541048959 
om  —  32.16174824838 
0(l)-3  .9921657339  28 

IER-  ERROR  FlAG-SEC  ASO  COMPUTER  SUBPROGRAM  LIBRARY  GUIOE 
IER-8 


COMPERISON  0*  ACTUAL  DATA  TO  XNTERPdLATEO  DATA 

A/M  -  CRA:k  LENGTH  RELATIVE  TO  WIDTH 
K-EYP  AO T UAL  OATA 
K-CAL  INTERPOLATED  OATA 
E  RELATIVE  ERROR  IN  X 


A/N 

K-CXP 

k-cal 

E 

.101% 

.9073 

.9467 

4.  34 

.1566 

1.0575 

.9776 

-7.56 

•  206% 

1.0572 

1.0545 

-.26 

.23/8 

1.1137 

1.1060 

-.69 

.2646 

1.1007 

1.1705 

6.35 

.2666 

1.1662 

1.1746 

-.98 

•  3096 

1.2366 

1.2590 

1.81 

.3460 

1.3490 

1.3244 

-1.82 

.3976 

1.3604 

1.4086 

2.04 

.6360 

1.4388 

1.4696 

2.14 

.6760 

1 . r  64  6 

1.5424 

-1.42 

•  4693 

1.6551 

1  .5710 

-5.08 

•  5760 

1.6481 

1.8774 

1.59 

.6467 

2.3460 

2.3688 

.97 

.7517 

3.9513 

3.9413 

-.25 
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COHPERISON  O'  ACTUAL  OATA  TO  INTERPOLATED  DATA 

A/W  -  crack  LENGTH  RELATIVE  TO  WIOTH 
ETA-EXP  A' T UA L  OATA 
ETA-CAL  INTERPOLATED  OATA 
E  RELATIVE  ERROR  IN  * 


A/M 

ETA-EXP 

ETA-CAL 

E 

1034 

1.9440 

1.9500 

.  31 

1556 

1.6643 

1.6491 

-.91 

206  4 

1.6855 

1.6334 

-3.09 

2  326 

1.6163 

1.6879 

4.43 

,26*6 

1.6950 

1.7933 

5.60 

,2666 

1.6796 

1.6011 

-4.16 

,3096 

1.9961 

1.9923 

-.  19 

,3460 

2.  1660 

2.1823 

.75 

,397# 

2.6055 

2.4945 

-4.76 

,4360 

2.7493 

2.7696 

.74 

,4760 

3.0632 

3.1251 

2.02 

4593 

3.2666 

3.2646 

-.66 

,5760 

4,5933 

4.6431 

1.06 

,6467 

6.6744 

6.5725 

-.76 

,7517 

12.0981 

12  1062 

.06 

THE  FOLLOIMG  “PINT  OUT  IS  FROH  THE  NUHERICAL  INTEGRATION 
SUBROUTINE (0:iDRE)  AND  INOICATS  THAT  SINGULARITIES  WHERE 
HANDLED 

I ER ’ERROR  FLAG, IF  65  OR  66  -SINGULARITIES  WHERE  SUCCESSFULLY  HANOLfO»ANY  ATMER  VALUE  CONSULT  INSL  MANUAL 
X1I)-THE  CRACK  LENGTH  FOR  WHICH  THE  ERROR  FIAG(IER)  IMPLIES 

GKHSTANCE  FROH  CRACK  TIP  RELATIVE  TO  WIDTH  (T/W)  FOR  WHf  C  H  THE  SINGULARITY  WAS  D£  f ACTED 


•••  WARNING  WITH  FIX  ERROR  (I £R  -  651  FROH  XMSL  ROUTINE  OCADRE 

X ( 1) «.  I 
Q-t.E-ti 
I E  R  *  65 

•••  WARNING  WITH  Fix  €PBj  R  (IfR  *  65)  FROM  I HSL  ROUTINE  OCADRE 

x  m».  i 
0-1  .£-12 
IER-65 


•••  WARNING  WITH  FIX  ERROR  (IEP 

X <Z) *. 15 

Q-l.E-10 

IER-65 

•••  WARNING  WITH  FIX  ERROR  HER 
x  <2) «.  15 
Q" 1 • E- 11 
IER-65 

•••  WARNING  WITH  Fix  ERROR  UER 

X (2)  «, 15 

a-i.E-12 

IER-65 


«  65)  FROM  IHSL  ROUTINE  OCADRE 

■  65)  FROM  IHSL  ROUTINE  OCAORE 

■  65)  FROM  INSL  ROUTINE  OCAORE 
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•••  WARNING  WITH  FIX 

X(3)«.2 

Q*  1  •  £-  10 

IER-65 

WARNING  Ml TM  FIX 

xm»,  z 

Q-l.E-U 

UR-65 

•**  warning  with  fix 

X(3)».2 

Q-l.E-12 

IER-65 


•••  WARNING  WITH  FIX 

X  (<*)-.  *5 

Q-l.E-10 

IER-65 

•••  WARNING  MITH  fix 
X  UJ-.25 

o-i.e-  u 

IER-65 

WARNING  with  fix 
X(WI  »,  25 
Q-i.e-12 
IER-65 


•••  warning  with  fix 

X 151 -« J 

O-t.E-9 

IER-65 

***  warning  WITH  fix 

X (5)*. 3 

Q-l.E-10 

IER-65 

•••  WARNING  with  fix 

X (5) -4 3 

Q-l.E-lt 

IER-65 

•••  WARNING  WITH  Fix 

X (51 -♦ 3 

Q-l.E-12 

IER-65 


warning  with  fix 

X  ( 6  )  « •  35 
Q-l.E-9 
IER-65 

***  WARNING  WITH  FIX 

X  (6)  -.  35 

Q-l.E-10 

IER-65 

warning  with  fix 

X  ( 6  )  -  .  35 
O-l.E-11 
IER-65 

•••  WARNING  with  fix 

X  (6)  -.  35 

Q-l.E-12 

IER-65 


ERROR  (IER  ■  65)  FROM 

ERROR  (ICR  •  65)  FRON 

ERROR  UER  -  65)  FRON 

ERROR  (ICR  •  65)  FROM 

ERROR  (IER  «  65)  FROH 

ERROR  (1 ER  •  65)  FROM 

ERROR  (UR  »  65)  FROM 

ERROR  (IER  -  65)  FROH 

ERROR  (UR  -  65)  FROH 

ERROR  UER  •  65)  FRON 

ERROR  UER  -  65)  FROH 

ERROR  (IER  *  65)  FRON 

ERROR  UER  -  65)  FROH 

ERROR  UER  -  65)  FROH 


IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  DCADRE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  DCADRE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  OCADRE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  OCADRE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  DCADRE 

IHSL  ROUTINE  OCAORE 

IHSL  ROUTINE  DCADRE 
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•**  HArNINU  KitN  /i*  tKhUR 

X  (7)  *.6 


G*  1 .  E-  9 
IER-65 

WARNING  KITH  FIX  ERROR 
X  (7  >«.<» 

0*1  •  E- 10 
IER-65 

•  ••  warning  kith  fix  error 

X  €  7  I  *.6 
Q*l.E-ll 
I  ER*65 

•••  WARNING  KITH  fix  error 

X  (71  ..6 

Q-i.f-i? 

XER*65 


UtN  *  OOt  TRU>1  WMliHL  UI44u»I. 

(IFR  •  691  FROM  IMSL  ROUTINE  0C4DRE 

(IfR  •  651  FROM  IMSL  ROUTINE  DCAORE 

(IER  -  65)  FROM  IMSL  ROUTINE  OCAORE 


•••  WARNING  WITH  Fix  ERROR  (HR  »  65)  FROM  IMSL  ROUTINE  OCAORE 

XU)«.65 
Q*  1  •  E-  9 
I E  R  *  65 

**•  WARNING  WITH  FIX  ERROR  (IER  -  65)  FRON  IMSL  ROUTINE  OCAORE 

X  (6  >  *.65 

o-i.F-ta 

IER-65 

•••  WARNING  WITH  FIX  ERROR  HER  •  65)  FROM  IMSL  ROUTINE  0C4DRE 

X  (5) *.65 
Q-t.l-ll 
IER- 65 


•••  WARNING  with  fix  error 
x (91 ».S 
Q* 1 •  E-9 
If R -65 

•••  WARNING  WITH  Fix  ERROR 
X  (91  *.5 
Q*t,E-10 
IER  »65 

•••  WARNING  WITH  FIX  ERROR 
X  (9)  *.5 
Q*l  .E-tl 
IER* 65 


(IER  -  65)  FROM  IMSL  ROUTINE  OCAORE 

(IER  •  65)  FROM  IMSL  ROUTINE  OCAORE 

(IER  •  65)  FROM  IMSL  ROUTINE  DCAORE 


WARNING  WITH  FIX  ERROR  (IER  •  65)  FROM  IMSL  ROUTINE  OCAORE 

X  ( 10)«.55 
0*1. e-9 
IER  *65 

•  M  WARNING  WITH  FIX  ERROR  (IER  •  65)  FROM  I  H^L  ROUTINE  OCAORE 

X  <10>-.55 
<JM  ,E-  10 
I E  R  *  65 

•**  WARNING  WtTH  FIX  ERROR  IIEP  ■  65)  FROM  IMSL  ROUTINE  OCAORE 

X  (10)-. 55 
QM  .f-li 
ICR *65 


WARNING  with  fix  ERROR  (IER  *  65)  FROM  IM^L  ROUTINE  DCAORE 

X  (U>*<6 
0*1.  E-9 
IER«65 

•••  WARNING  WITH  Fix  ERROR  (IER  »  65)  FROM  IK*L  ROUTINE  0C4CRE 

X  (ll)*«6 
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IER-65 

•••  warning  with  FIX  ERROR  IIER  ■  *5)  FROM  IHSL  ROUTINE  OCAORE 

xui)*«6 

Q-l.C-U 

IER-65 


•••  WARNING  WITH  FIX  ERROR  HER  ■  451  FROM  IHSL  ROUTINE  OCADRE 

X  II2)*.6S 

Q«t.e-a 

I ER-65 

WARNING  WITH  FIX  ERROR  HER  ■  451  FRON  I MSL  ROUTINE  OCAORE 

X<12>*.65 
0*1. E-5 
IER-65 

•••  WARNING  WITH  Fix  ERROR  (ICR  •  45)  FRON  IHSL  ROUTINE  OCADRE 

X<12>».65 

Q-l •£- 10 

IER*65 

•••  WARNING  WITH  FIX  ERROR  OCR  ■  45)  FRON  IHSL  ROUTINE  OCAORE 

X  < 1 2)* .45 

Q-l.E-U 

ICR-65 


•  WARNING  WITH  FIX  ERROR 
131*. 7 

1  .E-7 
9-66 

•  WARNING  WITH  FIX  ERROR 
13)*.  7 

t.e-a 

9-66 

•  WARNING  WITH  FIX  ERROR 
13I-.7 

l.E-9 

9-66 

•  WARNING  WITH  FIX  ERROR 

iJ)-.r 

l.C-IO 

9-66 

•  WARNING  WITH  FtX  ERROR 
1 3) * .7 

1  .C-ll 

9-66 


WARNING  with  FIX  ERROR 
X<l<*>-.75 
Q-l.E-8 
ICR-65 

•••  WARNING  WITH  FIX  ERROR 

XUAI-.75 

0*1.  E-9 

ICR-6? 

•••  WARNING  WITH  FIX  ERROR 
X ( l A )- .75 
Q-i  .€*10 
I E  R  *  65 

•••  WARNING  with  FIX  ERROR 

X I1A>* .75 

Q-l.C-U 

ICR-6? 


•  66)  FRON  IHSL  ROUTINE  OCADRE 


■  66)  FRON  IHSL  ROUTINE  OCAORE 


•  66)  FROH  IHSL  ROUTINE  OCADRE 


*  66)  FRON  IHSL  ROUTINE  OCADRE 


•  66)  FRON  IHSL  ROUTINE  DCAORE 


■  65)  FRON  IHSL  ROUTINE  OCADRE 


•  65)  FRON  IHSL  ROUTINE  OCADRE 


•  65)  CRON  IHSL  ROUTINE  OCADRE 


•  65)  FRON  I -SL  ROUTINE  OCADRE 
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MU'  |  ***  LI  «MH 

X (15I-.6 

Q-l.E-7 

IER-66 

•••  WARNING  WITH  FIX  ERROR 

0-1. £-6 
IER-66 

•••  WARNING  WITH  FIX  ERROR 

X(15>««4 

0*1.  E-9 

ICR-66 

•••  WARNING  WITH  FIX  ERROR 

X(15)-.6 

0- 1  *  E-  10 

UR-66 

WARNING  WItM  FIX  ERROR 
X<15)-«6 
Q  » 1  •  E- 1 1 
ICR-66 


(ICR  ■  661  FROH  IMSL  ROUTINE  OCAORE 

(UR  ■  66)  FROH  IHFL  ROUTINE  OCAORE 

(ICR  •  66)  FROH  1 HSL  ROUTINE  OCAORE 

UER  •  66)  FROH  I HSL  ROUTINE  OCAORE 


RESULTS  FOl  K  UNIFORM  TENSION  USING  THE  WEIGHT  FUNCTION  HETH00( K-NON  DIHENSIONAL 


A/W  -  CRACK  LENGTH  RELATIVE  TO  WIOTH 

K-EXP  •  VALUES  CALCULATED  THROUGH  THE  WEIGH  FUNCTION 

K-CAL  -  SJPPLIEO  COHPERI  SON  VALUES 

E  -  RELATIVE  ERROR  IN  X 


A/W 

K-EXP 

K-CAL 

E 

iooa 

1*2705 

1.1957 

6.  26 

19  00 

1.2697 

1.2662 

.  12 

2000 

1  *  WO  6 1 

1.3667 

2.74 

>2900 

1.5540 

1.4961 

4.01 

3000 

(.6059 

1.6551 

1.63 

>3900 

1.6170 

1.6965 

-2.  13 

6000 

1.9967 

2.1060 

-5.26 

>6900 

2.2917 

2.4  241 

-5.46 

,9000 

2.7719 

2.6266 

-1.93 

9900 

3.W792 

3. 3406 

3.76 

>6  000 

4. 36  36 

4.0632 

7.92 

6900 

5. 3196 

4,9993 

6.33 

>7000 

6.1632 

6.3759 

rS.02 

7900 

6.6662 

6.4009 

-19.04 

>6111 

7.3606 

11.9926 

-36.79 
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APPENDIX  E 

CRACK  SHAPE--CONICAL  APPROXIMATION 

Following  Orange  (Reference  11),  the  shape  of  a  crack  can  be 
approximated  by  a  conic  section  equation.  Requiring  that  the  crack 
opening  at  the  crack  mouth  will  match  the  actual  crack  mouth  opening 
n0,  one  can  write  a  general  form  for  the  conic  section;  equation; 


Figure  E-l.  Crack  Opening  Conical  Approximation 

where  t  is  the  distance  from  the  crack  tip,  a  is  the  crack  length,  and  m 
is  the  conic  section  coefficient  to  be  determined.  The  type  of  conical 
section  depends  on  m  if 

in  <  0  -  ellipse 

m  =  0  -  parabola 

0<m<°°  -  hyperbola 

m  =  ®  -  a  pair  of  straight  lines 
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In  order  to  determine  the  unknown  conic  section  coefficient  m, 
match  the  radius  of  curvature  at  the  crack  tip  from  the  Equation  E-l  to 
the  curvature  resulting  from  the  crack  tip  displacement  field  as  expressed 
in  Equation  16. 


Using  the  following  definition  for  the  curvature  of  a  two-dimensional 
curve  defined  by 


y  =  f(x) 


(E-2) 


the  curvature  K  is  given  by 


d2x 

dy2 


and  the  radius  of  curvature  is 

"■[11 

provided  K  ^  0 


If  one  uses  the  near  crack  tip  displacement  field  from  Equation  16, 


3t  _  tth! 

an  4KX 


n 


(E-5) 

(E-6) 


and  at  the  crack  tip 


AFWAL- TR-80-400 1 


From  Equation  E-6 


Also  from  Equation  E-6 


9 1 

an 


n=o 


o 


9  2 1  _ 
9n2 


(E-7) 


(E-8) 


Substituting  Equations  E-7,  E-8  into  Equation  E-3  and  then  using  Equation 
E-4,  the  radius  of  curvature  at  the  crack  tip  from  the  crack  tip  displace¬ 
ment  field  is 


R 


(E-9) 


Now  find  the  radius  of  curvature  from  the  conic  section  Equation  E-l 


n2  = 


2n. 


mn 


(2+m) a 


-  t 


( 2+m) a2 


(E-10) 


Differentiate  twice  with  respect  to  n 


2n 


2n  = 


( 2+m) a 


a  2mn 
9 1  o 

9n 


(2+m) a2 


9t 

3n 


2n 


3  2  t 


2  = 


( 2  +m )  a  3  n 4 


+  r  (hy+  t  at ] 

(2+m)  a2  LX50/  3n2J 


(E-l 1 ) 


(E- 12) 


From  Equations  E-ll  and  E - 1 2  at  the  crack  tip 

t=0,  n=0 
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The  value 


and 


5 2 1 
8h2 


( 2+m)  a 


(E-13) 


If  Equations  E-3  and  E-4  are  used,  the  radius  of  curvature  at  the  tip  of 
the  conic  section  is 


(m+2) a 


(E-14) 


Now,  matching  the  radius  of  curvature  at  the  crack  tip  of  the  conic  section 
(Equation  E-14)  to  .he  one  that  is  derived  from  the  crack  tip  displace¬ 
ment  field  Equation  E-9,  one  can  determine  the  conic  section  coefficient 


fej-  • 


Thus,  following  Orange  and  substituting  the  stress  intensity  factor  Kj , 
with  the  nondimensional  stress  intensity  factor  V  defined  by  the  relation 

Kj  =  ya/a 


The  value  of  m  becomes; 


2 


( E  - 1 6 ) 


where 


Y 


K 

o/a 


( E  - 1 7 ) 
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